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SYNOPSIS 


V j j} 


Thesis entitled, ' Some Stxjdies of Dynamical Symmetry Breaking 
in Gauge Models’ ty Janardan Prasad Singh in partial fulfilment 
of the requirements of the Ph.D. degree of the Department of 
Physics, Indian Institute of Technology, Kanpur. 

The idea of generating the masses of fundamental 
particles by their self -interactions has drawn a great deal 
of attention after the development of Weinberg-Salam model of 
electro— weak interactions. One way to do this is via the 
“Higgs mechanism". But this procedure is neither theoreti- 
cally appealing nor is backed by any ej<perimental ‘evidence. 
Moreover, it creates theoretical obstacles; asymptotic free- 
dom rrey be easily spoiled, there may be difficulties with 
quantum gravity, and there is the "fine tuning" problem. 
Alternative to this is dynamical symmetry breaking (DSB). 

Here, the asymmetry of the vacuum manifests itself through 
the non- vanishing vacuum expectation value of a composite 
operator and the scalar bosons occur as bound states. 

The mass scale in the originally massless theory is set 
essentially by the process of renorroaliaation. However, 
whatever one gains by adopting DSB, one pays for by loosing 
control on the theory (its solutions) and general ability to 
perform even approximate calculations. In view of this, it 
seems imperative to do calculations even 4,n simple models 
with more rigour. 

In the present thesis, we re-examine the Abelian 
gauge model proposed by Jackiw and Johnson, where a massless 



fermion field has a vector and an axial-vector gauge inter- 
actions. It has been claimed that DSB is realized in this 
model. The ratio of the axial-vector- boson mass to the 

fermion mass ('^) has been calculated by different authors 
m 

employing different approximations. We have investigated this 
model using the ’’gauge technique" which is suitable for low- 
-energy calculations such as the present one. In the limit 
where the ratio of the axial to the vector coupling constants 
becomes small ^ or, consistently, in the limit where ^ becomes 
small, we have derived an approximate solution for the fermion 
spectral function. This, in contrast to the findings of 
several authors, gives an extremely small value for ^ . 

Non-Abelian gauge theories, being asymptotically 
free, are better candidates to look for DSB. The strength of 
the interaction at long ranges in such theories can become 
very large. Given an attractive channel such forces might 
produce bound states which could act as dynamical Goldstone 
particles. Among the realistic cases where DSB has been 
applied, most successfiol one is the chiral symmetry breaking 
(XSB) in QCD. XSB solutions for the quark propagator have 
been used to calculate various quantities in low-energy 
hadron physics, such as pion decay constant, pion electro- 
magnetic form factor, amplitude for ^ yy ^ etc., and also 
to resolve u(l) problem in QCD. The results obtained are 
encoxiraging. 

We further apply these XSB solutions to calculate 
anomalous magnetic moments (a.m.m, 's) of up- and down-quarks 



and their mass difference. These are the modern version of 
old problems of a.m.m.'s of proton ansS neutron and their nass 
difference if one seeks explanation in terms of the consti- 
tuent quarks. 

To calculate a.m.m. of light quarks, we have used 
XSB solutions of quark propagator of both kinds: solutions 
based on infrared and asymptotic properties of QCD. Calcxila- 
tions have been done in one gluon exchange approximation using 
both "singular” and "nonsingular" forms of gluon propagator. 

We find that "singular" form of gluon propagator in axial 
gauge, and the XSB solution obtained by using the above 
propagator and longitudinal part of quark-gluon vertex form 
a good combination. The quark a.m.m.'s obtained in this way, 
along with effective qmrk mass, have been used to calculate 
proton- and neutron-magnetic moments using regxiLar SU(6) 
state vectors for baryons. We find a good agreement with 
experimental results. With certain assumptions, it has also 
been used to calcxiLate both spin dependent potential between 
a heavy quark and a light antiquark, and the differential 
cross-section for qq qq. In these cases, no significant 
effect of the a.m.m, term over the charge term has been found. 

In the gap equation for dynamical quark mass 
we introduce C3ED. Using asymptotic solution for 2^ for the 
unperturbed case, we find that to the lowest order in coup- 
ling constants the wrong sign of 2^ ” persists for all 
reasonable values of momenta. 

Our work strengthens the belief that mass generation 
is unlikely in case of completely massless Abelian gauge 



mc3dels. Our calculation of a.m.m. of light quarks and some 
of its phenomenological applications emphasize the role of 
dynamical quark mass for hadron physics. 
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Chapter 1 
INTRODUCTION 

"I 

With the recent discoveries of w“ and Z°, all aspects 
of the standard SU^(3 )xSUj^(2)xU(1 ) model seem to be well 
established with the exception of the Higgs sector, Hven if 
the standard model is integrated in a grand unified model 
such as SU(5)j one is left with an smbarrassing number of 
free parameters the self couplings and Yuk;aw.a couplings of 
the Higgs bosons with which to fit the data. In all 
modern fundamental theories, Higgs fields are introduced 
sorely to effect the breakdown of symmetry (including global 
chiral symmetry). There is no other compelling theoretical 
reason for scalar fields (save supersymmetry). Moreover, 
their pre-aence creates theoretical obstacles; They can easily 
spoil the asyomptotic freedom that was at the heart of gauge 
coupling unification. The reQ^uirement of the perturbation 
theory for the weak interaction that there must be at 
least one physical Higgs boson with = 0 (m^), and in 
any case i 0 (1) TeV, Is difficult to meett An elementary 
scalar particle can acquire large masses = 0 a-nd 

OCm^) where nip is the Planck mass Io''^ ^nd is the 
grand unification mass scale 10^ ^ GeV) from its propagation 
through short distance structure of space-time and from its 
interaction with heavy Higgses (with masses m^), while 
the radiative corrections will yield 6m| = 0(a^)x0(A^) which 
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is difficult to adjust with cut-off « mp or 

Alternatively, one may look for the source of spontaneous 
symmetry breaking (SSB) in the dynamics of the gauge theories 
themselves. If the asymmetry of the vacuum manifests itself 
through the non-vanishing vacuum expectation value of a 
composite operator rather than an elementary scalar field 
present in the lagrangian, then this is referred to as 
dynamical symmetry breaking (DSB), Here the scalar bosons 
(in particular, Goldstone bosons) occur as bound states. 

The massless Goldstone bosons may appear as physical particles 
on account of originally massless fermions acquiring mass 
(chiral symmetry breaking), or they may combine with the 
gauge bosons to make them massive (gauge symmetry breaking) 
exactly like the Higgs case. 

naturally, in such theories there will be fewer free 
parameters. In fact, if one supposes that the world is 
described by a grand unified dynamically broken gauge theory, 
then the single gauge coupling can be traded for the mass 
scale using dimensional transmutation. The mass scale in 
the originally completely massless theory is set essentially 
by the process of renormalization, necessary to remove the 
divergences present in any renormalizable theory. In such 
a world all mass ratios and the fine-structure constants are 
in principle calculable. The price for such an economical 
choice may be that in practice any calculation worth a 
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practical significance will be extremely difficult. 

The question of whether a given symmetry is broken or 
not is a dsmamical one. In theories with Higgs scalars one 
can break any symmetry by manipulating the Higgs lagrangian. 
DSB provides a self-consistent way to realize SSB, 

In case of DSB, there is no mass renormalization 
because there is no mass to be renormalized. The softness of 
dynamical mass generation is to be distinguished from the 
generation of mass from Higgs fields acquiring a vacuum 
expectation value. In the latter case mass renormalization 
is required and the fermion mass is not, in general, 

3 

computable . There is the further suggestion, at least in 
some approaches to DSB (e,g. in technicolor theories) that a 
new region 100 GeV of high-energy physics with heavy 
quarks and leptons must appear. 

In DSB, the scale of symmetry breaking is the scale 
at which the dynamics beccme strong. There are no dimensionful 
couplings, as there are in scalar theories, to set the scale 
of the dynamics. The dimensionless coupling constant of a 
non-Abelian gauge theory varies with the logarithm of the 
energy and it may take a huge change of energy before the 
coupling becomes large. For these reasons DSB is much more 
attractive in grand unifiod theories than the usual symmetry 

A 

breaking by scalar fields . Besides this, DSB figures 

prominently in the recent surge of interest in supersymmetric 

5 

theories , 


4 


The idea of DSB was borrowed originally from the 

g 

theory of superconductivity , where a transverse photon 

acquires a mass My by "eating up" electron -hole bound 

state such tho.t the free magnetic field equation v B = 0 

7 

becomes 

2 2 X 

V B = ^Jy B or ~ e (1.1) 

which explains Meissner damping qualitatively. 

There is a simple, intuitive , necessary condition for 
the occurrence of DSB in a given channel: The forces in 
that channel must be sufficiently attractive to bind a 
zero-mass Goldstone boson. In particular, in asymptotically 
free theories the strength of the interaction at' long ranges, 
or at small momenta, can become very large, irrespective of 
the value of the "physical" couplings. Given an attractive 
channel such forces might inevitably produce bound states 
which could act as dynamical Goldstone bosons^. 

The existence of any bound state, and in particular 
of a zero -mass bound state, can be examined by a Bethe- 
Salpeter (BS) equation. The existence of e solution of 
the BS equation is equivalent to the existence of a ay mm ia t.-ry 
breaking solution to the Schwinger-Dyson (SD) equation of 

9 

the theory , Usually the existence of the symmetry breaking 
solutions is not shown to occur necessarily, but rather 
that the integral equationsof the theory are consistent with 
such solutions. One assumes here that the nonperturbative 



masses depend analytically on the renormalized coupling 

constants, but that conTentional perturbation theory will 

not reveal this structure. Another possibility is that 

the mass depends nonanalytically on the renormalized coupling 
10 11 

constants ’ . The energy gap in a superconductor has 

similar nonanalytic behavior in the electron-electron 
1 2 

coupling 

The above method in linearised approximation yields 
a homogeneous equation which admits trivial solution (no 
symmetry breaking) as well as the symmetry breaking solution. 
Hence, one is unable to determine whether SSB does in fact occur 

In another approach, one calculates the effective 

potential for composite operators, and that enables one to 

study the stability of the syst®Q in a way that closely 

parallels the traditional treatment with fundamental scalars 

1 3 

a.nd to develop a stability criterion for candidate vacua , 

. 1 4 

But as recently shown by Miller , an effective potential 

evaluated perturbatively can possess false minima, artefacts 

of the approximation itself. In the present thesis, this 

approach will not be adopted. 

Of course, whatever one gains by adopting DSB one 
pays for by loosing one's control, on the theory (its 
solutions) and, short of acquiring a considerable amount of 
insight into nonperturbative effects, one's general ability 
to perform even approximate calculations. In view of this 



fact, it is important to do calculations even in simple 

models with more rigour. The experience gained from such 

calculations can be used later for more complicated and 

realistic cases. In Chapter 2, we shall introduce an early 

1 5 

Abelian ga-uge model proposed by Jackiw and Johnson , and 
then investigate it more rigorously using the so called 
"gauge technique”. In Chapter 3, we review chiral symmetry 
breaking (XSB) in QCD by deriving the gap equation and 
finding out its asymptotic solution, XSB solution has also 
been founi out in infrared region through the SB equation 
using "singular” foim of the gluon propagator. We illustrate 
the use of dynamical quark mass by showing how pion decay 
constant is calculated . In Chapter 4, we use XSB 
solutions obtained by various authors to calculate anomalous 
magnetic moments (a.ra, m.’s) of light quarks. We use solutions 
of both kinds ; solutions based on infrared and asymptotic 
of QCl>. We also work out some phenomenological 
implications of the presence of a.m.m. In Chapter 5, we 
calculate u-d quark mass difference and then comment on 
possible reasons for getting wrong sign of the mass 
difference, Finally, Chapter 6 contains some concluding 
remarks. Some of the mathematical details have been put in 
Appendices, 
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Chapter 2 
A MODSl SXAlvIPiS 


Introducti nn 

I^SB in particle physics was initiated by the 
Nainbu-Jona-lasinio model^ based on four-Permi interaction 
where it was shown that a theory may possess a XSB bound - 
state Goldstone boson. Next Schwinger^ made the crucial 
observation that it is possible to have a pole in the polari- 
zation tensor of a vector gauge boson at zero momentum 
transfer. This pole leads to SSB by giving a mass to the 
otherwise massless gauge particle, a process known as the 
Schwinger mechanism. It turns out that the Schwinger mechanism 
is unique for 1+1 dimensional and it occurs due to 

topological reasons^. Another eirgument why the Schwinger 
mechanism does not work for, say, 3+1 dimensional QED is 
due to charge-conjugation and I-invariance'' If charge- 
conjugation is unbroken, a vector field with odd charge- 
conjugation can not absorb a Goldstone boson which is a 
bound-state of fermion-antifermion pair in the scalar state 
having even charge-conjugation. A similar argument works 
for I-invariance also. By contrast, an axial -vector gauge 
field is C even and T-rodd as is the bound pseudoscalar 
state of fermion-antifermion pair. So axial -vector and 
mixed parity gauge fields can acquire mass in this way. 
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But a pure axial -vet; tea? exchange is repulsive and can not 
fo 2 m a board -state. On the other hand if the theory contains 
both vector as well as axial -vector gauge particles and 
the vector coupling is stronger than the axial coupling so 
that the net attractive force can form a bound-state Goldstone 
boson, then one can hope that the Schwinger mechanism will 

1 5 

work for the axial -vector gauge boson. The Jackiw -Johns on 
(JJ) model implements this idea in an Abelian gauge model. 

It has been claimed that in this model both chiral symmetry 
as well as axial -vector gauge symmetry get broken , 

However, different authors have employed different approxi- 
mations and consequently have obtained different results 
for the ratio of the axial-vector-boson mass to the fermion 
mass. This is one of the reasons why we took up this model. 
What we are going to describe in the rest of this Chapter ie 
contained in Reference [21] , 

2,2 Review of Barlier Work • 

The JJ model is described by a zero bare-mass 
lagrangian given by 

X = ? (i;^ -g / - g' -1/40^ 

-l/4(3|iBa;- (2.1) 

which is gauge invariant under the full 11^(1 )xUj^(1 ) set of 
transform ations, namely the ordinary phase transformations 



(x)/g 


^ (x) -♦ exp (i A. (x) ) ^' (x ) ,, Ajj^-* 
as well as the chiral set 

'l(x) - exp(iX' (x) Y 5) 'i^(x) , B,, Ai (x)/g' . 

In searching for the phase which breaks chiral 

symmetry, namely <'if> 0, one assumes that the fermion 

acquires a non-zero (renormalized) mass and correspondingly 

that the inverse fermion propagator S (p) develops a 

non-?( component. This will imply, via the axial Ward- 

Takahashi (WT) identity, that the axial -vector vertex 

has a pole which is associated with the induced G-oldstone 

excitation. As a result the axial self-energy develops a 

1 5 

pole in its momentum variable. This in turn implies a 
non-zero axial -vector-meson mass |j, , Therefore, in such a 
chiral model, one can anticipate a relation between M'/m 
and the axial -meson coupling constant g’, where m is the 
dynamically generated fermion mass, Furthermore, it is also 
clear that \x should tend to zero as g' -» 0 (that is, in 
the free-field-theory limit). 

As emphasized by Delbourgo and Keck”'^ (DK), the gauge 
22 

technique may be a suitable tool for investigating such 
probleans since the gauge identities must be appropriately 
respected whether or not USB ensues. Secondly, USB is a 
low-energy phenomenon which occurs due to bound-state formation, 
and it is precisely this ^region of energy where the gauge 



23 24 

technique works very well ’ . Below, we shall briefly 
review the work of DK and fix up our notation. 

We hare vector as well as axial WT identities for 
the renormalized Green functions 

k*^S(p).ria (p,p-lc) S(p-k) = S(p-k) - S(p) (2.2a) 

k^S(p)r^5(p,p-k) S(p-k) = Y5 S(p-k) + S(p)if^ . (2*2b) 

Hace S den®tes the complete fermion propagator, whose 
iinre normalized form is given by 

iS(p)p^ = f y i (x-y) 

= /dV 0lT(^'p(x)fa (y))J 0>.. 

(2.2o) 

p's stand for the irreducible vertex functions with 
coupling conetants factorized out. Their unrenormalized 
forms are defined by 

iS(p) rp.(p,p' ) i S(p' ) = / d'^x d'^'y <0 I T 

('i'(x) Jp(o) f(y)) 1 0 > (2.2d) 

and 

iS(p) ni5 (P>P' ) iS(p') = / d'^x d^y e“^P' OIT 

(f(x) T (y))fO>, ( 2 .. 2 e) 


where 



and 


1 1 


Jp.5(x) = Y5 ^ (^) (2.2g) 

Since Equation (2.2b) has been derived in a formal way on 
the assumption of axial current conservation, one has to 
ensure that quantum loop corrections do not invalidate it. 

This is possible only if regularization anomalies are absent 
or else cancel—in practice this may mean a doubling of 
fermion fields. We have further assumed that Equation (2,2b), 
valid for unrenormalized Green's functions, retains its form 
after a multiplicative renormalization. 

If we write the fermion propagator in terms of a 
spectral function as 

S(p) = (T , ?) ( 2 . 3 ) 

m p-W + ie'£(W) 

where P (W) is the fermion spectral function (Appendix A), 
then, up to a transverse correction, an appropriate solution 
of Equations (2,2a) and (2.2b) can be written as ^ 

’ dWP (W) 

S(p)rp(P,P-k) S(p-k) =/- — r Y — (2,4a) 

S(p) ru5(P»P"^k) S(p-k) = f 2"-^) • 


(2.4b) 
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The appearance of the pole term at the axial leg indicates 
the presence of the pseudoscalar Goldstone mode. This 
induces a dynamical singularity in the meson self -energy 

= ig' Z2/(dp) 'TrCSCp) ^^c(p,p-k)S(p-k)Y^^ Ytr)»^dp)= — 

{ 2 %)^ 

(2;5) 

as k -♦ 0 and this, in turn implies that the original zero- 
mass axial meson acquires a finite mass p . Here Z 2 is the 
fermion wave function renormalization constant. Substituting 
Equation (2.4b) in Equation (2,5) we have 


knkj, 

(k) + (-Y- 


S^y)8 


ig' 


//■ 


(dp) dW W^P (W) 
(p^-W^)[(p-k)^-W^] 


(2.6) 


where is the conventional QED self-energy for a coupling 
g' . Invoking conservation of axial current, this can be 
written as 

’^^50^5 (1^) ^ ""a (2.7) 

then clearly 

■* xVk^' as 0 (2,8) 

by virtue of the Goldstone pole in 


Here 
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P (dp) dW W^P (W) 

X = -8 i Zp // 2 . (2.9) 

(p - w ) 

Now it can be seen that the renormalized inverse propagator 
(Zj being the axial-meson wave function renormalization 
constant ) 

•*1 2 2 2 2 

)(k z^+k 7i^(k ))+a gauge fixing term 

( 2 . 10 ) 

2 

no longer has a zero at k =0. It is shifted elsewhere, 

2 2 

to the physical (mass) P which can be determined by invoking 

—1 2 2 2 2 

the renormalization condition D (k -F) near k =p . In 
the pole approximation 

2 ,2 2 2 

p, =g \ =“8ig* Zp 2 • (2,11a) 

(p - VT ) 

2 ‘^2 

After Wick rotation (p_-*- iP/? p = Z p ) and introducing 

0 4 i=1 ^ 

the dimensionless quantities 

0 )= W/m, r(w) = mP(W) , ^ = p^/m^, 

Equation (2.11a) can be written, after the angular integration, 
in the form 



(2.11b) 
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To determine r(a)) or p(W), one writes the Schwinger-Dyson 
(SD) equation for the fermion propagators 


-1 , P ^ 4^ 

Zg = S(p) ^-±g / (dk) S(p)rp. (p,p-k) S(p-k) Dj^ {^)yy 

-ig' f (dk) S(p) r^5(p,p-k)s (p-k) Tg (k)yyy^ , 

( 2 . 12 ) 

The complete gauge boson propagators in Sq,(2,l2) are 
replaced by the bare ones, which in the landau gauge can be 
written as 


\xv 


+ie) 



- »^A,B + 


^A = ^"b = <^2.13) 


Then using Equations (2,4a) and. (2.4b) in Equation (2,12) and 
.the fact that (Appendix A) 

= / P (W) dW 


one arrives at the following self-consistent equation for 
P(,W) (derived in Appendix A)s 


-£(W).P(W)(W-m) =/ Im (yw, W ) +gj(W,W')) 

(2.U) 


Where (see Appendix B): 


Im S»(W,W') = (w^-i^!^)©(W^-W'^) 

. 16 T^\r 


(2.14a) 



15 


and 


,2 o 2 2 (W^-W'^)^ 

ImST.(W,V/') = ((r+6WW>+W>^-2^ )0 + (0,^ -0J) 

"-T rs T r-^ “ 


32 ^^W 




1-i 


(2.14b) 


1 /2 

0 ,^ ={[w^"(W'+iA)^j[w^-(W'-!.i)^] } e&v’^-(W’ + i-i)^]. 


In view of the fact that the integral equation (2.14) is 
intractable, DK resorted to an assumed form of P (W): 

m p(W) = 8(W) ;:3- ;; [ | F( -S- dt ^ -6- 6^ ; -26 • 1 - ^) 


r (- 6 fd ') r (- 26 ) 


m 


+ (-1±S.) P(-6-6>, 1-6-fii; -26; 1-4)] 


6 - 6 ' 


m 


with 6 = 3g^/l67t^, = 3g’^/l6 , This guessed form of 

P(W) satisfies the following requirements: (i) it has the 
correct infrared and ultraviolet behavior and (ii) it reduces 
to the correct QED form in the limit g’ -* 0. Their guessed 
form of P (W) gives the following result: 


In a later paper Delbourgo and Kenny claim to have made 
a better guess for P (W) which we have exhibited later in 
this Chapter, But the result for [i/m comes out to be the 
same in both cases. 
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2,3 Approximate Solution of P (W) and Equation for the 
Mass Ratio 

Let us consider a limiting case when g' <<g. Then 

2 2 

DK's result for ^i/m is 0(g' /g ) whereas JJ obtained 

o ^ 2 

2/2 4 S' 

H /m £1 3 -2“^' ^ > 

g “ g' 

which, in the above limit, will give P-Zm to be 0(g'/g)» 

Thus we can attempt to obtain the value of P-Zm in that 
limit since we can hope to solve the integral equa-tion in 
that unit: we shall employ a solution obtained by 
(non-perturbative) gauge technique in the gauge coupling g 
and by a perturbation expansion in g' , We will look for a 
consistent solution with p.Zm<<1. Needless to say, the 
P (W) determined in this way should satisfy the requirements 
stated in S ection 2, 

In the limit stated above, one can attempt an expansion 

of P (W) in a power series in \i/vi t P (W) =:P°(W) +AP (W), 

Here P° (V/) is the fermion spectral function with p = 0 and 

AP (W) = (i-i/ m) P''^^ (W) )P^^^ (W) + 

It can be easily checked^ as will be shov.n later, that 

22 

P°(W) will have the same form as found in pure QBD with 

2 2 2 o 

the replacement g ^ g g' * When this P°(W) is used to 
calculate the axial r-meson mass, one gets a zero answer 
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(shown in Appendix C) as one should, since the Goldstone 

1 9 

mode does not materialise in QSD . It has been shown in 

References [23 J and [24 ] that the fermion (interacting with 

a massless gauge boson) spectral function determined from 

the longitudinal vertex is sufficient to determine the infrared 

properties of the thecryi hence the correction to p'^(W) 

coming from inclusion of the transverse vertex would not 

give a non-zero value for p, . Thus the gauge technique 

tells us that in all calculations where the axial -meson 

propagator with 0 has been used to calculate the fermion 

mass (m/^0) , and subsequently the axial -meson mass itself, 

one should get p,=0. Thus it is only A P (W) which can give 

a non-zero contribution, if any, to p . Our attempt will 

be confined to calculate caiLy the first term in the expansion 

of A P (W). Since the coefficient p/m of (VT) comes 

( 1 ) 

from an internal meson propagator, P' ' must have an overall 

2 ' ^ 2 

factor of g' as we shall show later, i,e,,AP (W) {|i/m)g’ 

in the lowest approximation. If this AP (W) is substituted in 
equation (2,1 lb), then one can see that the expression 
obtained for pi/m in this way will have a factor of g’^^ 

This is our main result. Below we shall give its derivation 
leaving some of the mathematical details fpr the Appendices, 

2 2 

We shall solve Equation (2.14) for g’ /g << 1, or, 
consistently, for \x/m «1, Our approximation consists of 
simplifying the expression for Im 2-g and then solving for 
P (W) perturbatively in the small parameters stated above. 



18 


Let U.S define 3g^/l6i'^ =6 ^ 3 g’^/l 67 F = 6 ’ , 6 - 6 * = 62 , 
and break r(w) intc even and odd parts as 

r(co) = r_^(w) + r^(c4 (2.15) 

where 

r_^(w) =-^{r(co) +r(-a))}, r_(a)) = ■|•{r( 6 )) -r(-w)}. (2.15» ) 

Ilow, in the infrared region, where (IWI - 
Equation (2.1 4 .) can be written as 

0 )^ e (oj) (w -1 )r (w) = - / ( 6 a)' ^ )+ 6 » 

oj “ 0)' 6 (41 /m 

e(a)^-a>'^). (2.16) 


The solution of Equation (2,16) with the last term dropped 

p 'Z OtZ 

is known (derived in Appendix C) 


r‘^'(a)) = 

.r (-2 S) 2 


(oT'' E(-6^ -6. -26; l-co ) 


+ E(- 6 , 1 - 6 ; - 26 ; 1 - 0 ) )). 


(2. IT) 


On writing 


r^ (w) 


r°' + A r^ (0)) 


(2.18) 


one finds (shown in Appendix D) that 


Ar; ((.)» - T 4 fii- 2 ^) ) r°'{to) 

— P. I — 


(2.19) 
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In a region where W is such that lJ'/(W-m) << 1 and 
remembering that in Equation (2,l4)jImS2 (W,W’ ) is multi- 
plied with the weight function P (VP ), which is concentrated 
mainly at W'=«i, Equation (2,14) can be written as (Appendix B) 


m^e( m) 
x(co' (o)^ 


(co-1) r(m) - ^ 

2 2 3 

-6)’^) + ) (^[^^^-(m’+ij/m)^ ] - ©( ctF-w ^ ) )J , 

6 ( 4 i /m 


(2.20) 


The solution of Equation (2,20) with u/m = 0 is easily 
found out to be 

e(a))(w2_i)-1-2.^2 -1 , , 2 

r°(ai) = c I'(-62, -62;-262;1-m ) 

r (- 262)2 2 

+E(-62, 1-625-262* 1-0?)), (2.21) 


where c is a constant which may depend on 6* , Again writing 


r(m) = r° ( co) +Ar(co) 

(2.22) 

one finds that (shown in Appendix D) 

A r^(u)=; - ^ S' (ii/m)(u2-i fl r°(u) 

00 ( 2 . 23 a) 

A r_(w) « - ^ 6 ' (^i/m)(co^-1 w) 

^4 “^.^ 00 ( 2 , 23 b) 
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where ^ is defined in such a way that and p,/m«(X-1) 

2 

(one can think of i-i/m = 0 [( X -1 ) ]) , 

We have one unknown constant c which could be 
determined by the condition that r(C(>) should be continuous 
at w= + (l4-p,/m). Unfortunately Ar(w) is not known for Iml 
just above 1+ii/m, However, it is clear from Equation (2,17) 
that 0 = 1+ 0(6'). let us compare our result for r_|_(m), 

2 5 

with that obtained by Delbourgo and Kenny (translated 
into our notation) 2 

r (a))~ £(01)0)“'' ( p(-6, -65 -26;1-«^) 

+ r(-26)2 

+ 26' E(1+ 6p, 1 + 6 5 1- 1-0?)). 

2 

In Table 1 we have tabulated both the results for 

2 

e(o))o) r,(o3) for different values of , It is to be noted 

2 

that in the calculation of (P-Zm) , it is only the second term 
in Delbourgo and Kenny's expression for r^(a)) which contri- 
butes, and in our case it is Ar^(u) and Ar_^(m), given 
by Equations (2.19) and (2,23 a) respectively, which contribute. 
These two expressions are different from each other. In 
our opinion our approximate form of P(W) is more reliable, 
in the region g' Zg <<1 and \i/m « 1 as it is based on 


* 


an explicit calculation. 



Table 1 

result for e ( to) wr,(w) with that obtained by Delbourgo and Kenny for 
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To calculate the mass ratio, as given by Equation 
(2,11b), we have to perform the following integration; 


^2 2s2 7 2 

(Q.4-Ct)J |Ci)l=1 (q+c^) 


CD 


+ // 


l oi I =1 + iym 


■2 


(q" 


First consider the integration over o) : 


-1 


^1+p/m (w) do) 

I ) — ?r"« — 


/•^x‘(Ci))dco //* /• - ,^1^/*-* j‘_! — 

/ — = C 5 

(q + ) -d-Hi/m) -1 (q + ) 


-(1+p/m) A 
+ [( J- + ^ 

1-ipi/ra 


-X 00 

) + ('^ + A ^ ] 

-CO 


A r_^(a)) do 


(2.24) 


The third term gives 

2 




00 


{ J" + /. ) d o 


-00 


K 


6’ (p/m)2 

(q^ .o^)'^ 9 r(-2 62) 


26, 


Ar^(co) 

1 


44 


1 

2 2 

(q +1 r 


2 -’-2*2 


• 4 <^2 -q^+i ■ 

2 ”1 2 n 

-I- terms of o [(^ ”1) J) 
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. M 1 , r^(1^a2r(-2a2)r(-1-262) 

9 r(-262) (1^+1 f (-62) 

-r(-262)r^(1+62) P(1+62, ) "(q^+l ) f (-262 ) 

xT^ (2 + 62) P (2 +6^, 2 +62? 25 - q^)). ( 2 . 25 ) 

In the first term of Equation (2.24) the contribution to the 
integral comes only from the limit u= + (1 + 11/m) . Ar_^(t)) 
in the range of the second integral is not known. However, 
if A r_^(a>) is smooth enough at co= + (1 tii/m) so that the 
integral is continuous at this point, then the whole 
contribution to Equation ( 2 , 24 ) comes from the second term 
of Equation (2,25)5 if the integral is not continuous at 
0)= + (1 +|i/m) then p /m will be logarithmically divergent. 
Assuming that the first possibility is realized, we get 

, X2 6.2 _ _i2 r(1-2^) 2 

" ■ 2'' V(2+262) 2(1+262) r"(1-«2) ^ 

( 2 . 26 ) 

where A is an ultraviolet cut-off. But in a truly consistent 


The behaviour of integrals involving, terms like r (m) in the 
infrared region has been discussed by Atkinson and Slim^T, 
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15 * 

theory of DSB l^/m should be finite . On the other hand, 
we observe that the logarithmically divergent term is 
smaller by a factor of 6 2 compared with the first term. 

It is at this order of magnitude in coupling constants that 
the contribution of the transverse correction to Ar_^(6)) 
becomes important if we assume that this contribution is 
non-zero in the (infrared) region of interest. If \i/m has to 
be finite then the divergent term in Equation (2,26) must 
cancel the term obtained from the transverse corrections to 
Ar_^(ci)) and we would get 

^/m = 0(S'^) . 


2,4 Conclusion 

Our result for the fermion spectral function agrees, 

25 

by and large, with that obtained by Belbourgo and Kenny 
in the infrared and ultraviolet regions, but in the inter- 
mediate region it disagrees with their result. This implies, 

in particular that we get the same asymptotic form of the 

•19 24 

fermion self-energy function as obtained by JJ ’ , But 


The consistency of our approach to this problem requires 
P'/m to be not only finite but small. 
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our result for 1^/m, in sharp contrast to the findings of 
earlier authors, is extremely small. 

The above result could be an indication that the 

dynamical gauge symmetry breaking hypothesis is questionable 

in the JJ model. Indeed using renormalization group arguments 
1 1 

G-ross has shown that in case of infrared free theories, 
as the interaction (one or more) responsible for the DSB is 
removed the masses increase. This, in particular, means that 
for an elementary particle mass one will not have a smooth 
transition to the free field theory limit. Gross regards 
this as an indication that infrared free theories can aot 
produce DSB, 

28 

stern , again using renormalization group arguments, 
has also shown that in completely massless Abelian gauge 
models, like the JJ model, dynamical mass generation is not 
possible, at least for weak couplings. However, Stern points 
out, unlike the case of the gauge boson, fermion mass gener- 
ation can occur at a nontrivial infrared -stable fixed point. 

We have worked in the landau gauge because in this 
gauge expressions for Im and Im are simple. In other 
gauges, it is difficult to find the solution for fermion 
spectral function even in QEP , like earlier authors, we 
have not shown the gauge independence of our result. 
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Since symmetry breaking in realistic physical models 
is expected to be an infrared -dominatei. phenomenon"^, we 
will now onwards concentrate on non-Abelian gauge models. 
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Chapter 3 

CHIRAL SYMviSTEY BRUIAKING IH QCD AND ITS APPLICATIONS 

3.1 Introduction 

It is becoming increasingly clear that in QCD there 
are two physically meaningful concepts of mass; the pertur- 
bative or current quark mass JOr- which occurs in the bare 
QCD Lagrangian (presumably due to weak interactions) and the 
nonperturbative or dynamically generated quark mass 
which dominates low energy hadron physics. To clarify it, 
consider the QCD lagrangian with m=Oj then it has 
SUj^(Nj,)xSUj^(Nj,) symmetry = number of flavors). The 
success of PCAC (partial conservation of the axial 
current) in conjunction with current algebra is an experi- 
mental fact. If we assume QCD describes strong interactions, 
this requires that the above symmetry is dynamically broken 

down to SUy (Nj,) : Q^ = 0, 0, a = -1,with 

2 * a 

Np,-1 attendant massless pseudoscalar bosons and Q =0, 

Q^I0> = 0, This Nambu-Goldstone (NG) realization of 
29 

chiral symmetry is the best place where a dynamically 
broken symmetry seems to be relevant for the real World. 

When chiral symmetry is explicitly broken by current quark 
mass (m ii^ 0) terms, the NG bosons become massive. 

To give an idea of why chiral symmetry is realized in 
the NG way, let us again use the analogy with superconductors. 
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The Cooper-pairing rearranges the surface of the Perm! sphere 
because electrons hare attractive interaction due to phonon 
exchange. In a similar way, the attraction in the scalar . 
qq channel is so great that it destroys the well-defined 
Dirac sea and produces, a finite gap (= effective mass) 
between negative and positive energy states. Another way 
to express this is to say that a quark condensate is 
developing! 


< 0 } ^(x) l(x) J 0 > 0, 


(3.1) 


also being analogous to the case in superconductors. To 

30 

give a numerical estimate, it is known that 


<0 I f T I 0> -1.6x10“^ GeV^' 

< 0 I (g 0 >~0.5 Ge?"^, 

where ^ is a quark field, and G^^j, is the gluon field 
strength. 


(3.2) 


In the literature, several *' causes" have been put 
forward for dynamical chiral symmetry breaking (XSB) in 
QCD. We mention a few of them. 

(i) Computer simulations in lattice gauge theory suggest 

31 

that the chiral symmetry breaks down spontaneously , 
and that the range of force responsible for XSB 

32 

is relatively short, independent of confinement , 
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Let us recall a crude but basic idea about how XSB 
occurs in gauge theories. Consider a bound state of 
a pair of massless qq. Because of t.ho uncertainty 
principle, the energy of the ground state (pion) will 
be given byS^ss p - g/r« p(1-g ) where p and r 
denote the relative momentum and coordinate, respectively. 
In a fully relativistic formulation, this relation may 
be replaced by « p -g yr » p (1-g^). When the gauge 
coupling constant exceeds order one, there will bo 
tachyon bound states, indicating instability of the 
vacuum. In order to cure this instability, the vacuum 
rearranges itself and gives mass to quark so as to 
keep the pion massless. Hence, in the ladder approxi- 
mation (or, more^ precisely, in the single gluon exchange 
approximation), XSB is expected to occur for coupling 
constants above a certain critical value. In fact a 
simple analysis of distance dependence of the coupling 

constant [a„(r) = t— ] shows that a (r)Rj 1 

b In ® 

at a distance of about Rw T^fm which is roughly the 
size of the hadrons. This sets the scale of non-TDertur— 
bative confining effects whereas other non-perturbative 
phenomena (such as ..XSB) start essentially already at 
distances of about 0,1 - 0,4 fra (depending on the 
particular quantity under consideration), where 
perturbation series still look good enough^^. 



(ii) Several authors”' have argued that XSB will 

occur, in anj^ theoiy of massless fermions interacting 

through non-Abelian confining gauge fields. Strong 

coupling lattice treatments precisely seem to demons- 

38 

trate at the same time XSB and confinement"^ , It 

has been proved, within the mean field approximation, 

that confinement and XSB are appearing or disappearing 
39 

together . 

1 6 

from a study of various models, Cornwall has 

illustrated how a theory of confinement automatically 

leads to XSB. In many such theories, the quark 

propagator is expected to be an entire function of 

the momentum, a characteristic criterion for confinement, 

A study of the SB equation for the quark propagator, in 

35 

a confining QCD, by Ball and Zachariasen reveals that, 
in addition to a chiral -symmetry preserving solution, 
there exists a broken symmetry solution which has no 
singularity, thus hinting at confiment. More recently, 

-7 C 

Acharya and ITarayana-Swamy , by an examination of 

WT identities and SB equations in QCB, were lead to 

the conclusion that there exists a connection between 

the NG realization of chiral symmetry and quark 

confinement. In fact, these authors have also 
40 

demonstrated that chiral symmetry remains unbroken 

due to single gluon exchange when the gluon propagator 

-2 2 

has the behavior B^j;~q as q **• 0 (representing 
nonconfining force), A similar connection between 



confinatnent and XSB has been shown in Reference [37]. 

(iii) Some authors hawe demonstrated that the instanton 

structure of QCD may break the chiral symmetry. The 
reason for this is, as proposed by Callan, Dashen and 
Gross , that instantons can generate attractiye forces 
between quark and antiquark* Since this topic is not 
directly relevant for our purpose, we shall not go 
into its details. It is to be pointed out that the 
quantum theory of interacting instantons turns out to 
be rather complicated, so the most interesting effects 
have not been so far reliably calculated from first 
principles, 

42 

However, in a model dependent approach , the 
instanton parameters, such as the instanton density 
n. in space-time and typical radius P. , have been 
fixed from data (mainly the gluon condensate)* Hrom 
these parameters, the quantitative estimates of < 
and m^^^ (quark effective mass) have been made. The 
numbers agree reasonably well with those obtained 
from other sources (e.g. equation (3*2) for > ). 
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3,2 G-ap Bquation and Its Solutions 

Here, we shall first derive a generalized gap equation 

2 

for the mass operator 2 (p ), resulting from the SD equation 

satisfied by the flavor-nonsinglet axial-vector vertex 

1 0 

function and the corresponding WT identity 

As emphasized byPagels"'^ and Cornvall^^, corresponding 

2 

to two concepts of mass, the mass operator 2 (p ) has two 

2 

distinct contributions: 2^ (p ,m) arises from the explicit 
breaking of chiral symmetry due to the current quark mass 
and 2^ = 0 when m = 0; the dynamical mass operator (P ) 
is independent of m and vanishes to any finite order in 
perturbation theory^^ Since our main object is the HG 

realization in the symmetric theory, we shall set m = 0 
from the beginning. In this case, 2 = 2^^ satisfies the 
gap equation. 

The qq Be the -S alp e ter scattering kernel K(p,q,k) 
(described in terms of truncated four fermion diagrams which 
can not be disconnected by cutting two fermion lines), 
which occurs in the gap equation, contains two distinct 
pieces: K = Kp + » where Kp admits the familiar skeleton 

expansion and denotes the nonperturbative contributions 
arising from instantons. K^jp vanishes to any finite order 
in perturbation theory. It has been suggest ed^^ that 

responsible for the HG realization of chiral symmetry 
in QCD, On the other hand, there also exist treatments of 
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the NG- realization of chiral symmetry based on the ladder- 
approximated kernel wherein the instanton contributions 
are totally ignored'^ 

The flavor nonsinglet, anomaly- free, renoimalized, 

Si 

proper axial-vector vertex satisfies the WT identity 

Z S' O' 

ir~X5 (PfP+'l) “ “ 2^ *- '' 5 "2” ^ (P+l) + ^ (p) Y 5 1 

(3.3) 


where we have defined the (unrenormalized) axial -vector 
vei'tex function by 

(iS(p)f"^^ (p,p+q)iS(p+q))^ = d^x dV 

5 a - 

x<0lT (s ¥^(x) ^^(y))}0> , (3.3a) 

r=1 / 

j,!, = i. 'ir ■ (3.3b) 

r = i 

Here the summation is over color indices, and JN.®' is the 

p 

flavor SU (N^) matrix with a = -1. Constants 

i *0/ 

Z2 and Zj^ renormalize (the unrenormalized) S and t~ 
respectively by multiplying them. The inverse of the" 
renormalized quark propagator has the form 

S"''(p) = A(p^) ^ - S^(P^) (3,4) 


and possesses exact flavor symmetry; 
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[s^^p), -^] = 0, a = 1 , 2 ,...,n 2 _ ^ ^ 

Assuming that SUj^(Np )xSt!^(Np ) spontaneously 'bi'ealcs down to 
SUy(N^), we have 

lim qVx5 (PiP+'l) S'"''(p)} ^ 

= (^)2 (p )y ^ (3.6) 

This corresponds to a hound-state pseudoscalar pole at 
2 ^ 

q = 0 in , the bound state transfoiming accoiding to 

the adjoint representation of SU(Nj,) and as a color singlet: 

(p,P+4)jp„ie = 4 *n<?(P.P+q)'<’5 -f . (3.7) 

Here_,{>>(p,p+q)Y 5 quark-quark-pion vertex function, 

and f^ the pion decay constant. We can combine Equations 
(3.6) and (3.7) to get the Goldberg er-Treitnan relation, 

z 

^tc^Cp^) = lim f7i;(?(p,P+q) = (^)2IL (p^), (3.8) 

q -*0 2 

The renormalized proper axial -vector vertex satisfies the 
SD equation 

r ^5 (p,p+<j) = Z^Y^Y 5 4- + [ 8 ( 1 :) r ^ 5 ( 11 :, k+q) 

X S(k+q) ] K(p,k,q) . g) 


When Equations (3,6) and (3-9) are combined together one 
gets the desired gap equation . 



^X) = - / (dk) [ ^ ^5 ^(p,k,0), 

(3,10) 


^or large momenta the perturbative contribution 
to the kernel Is dominant and we can Ignore the nonpertur- 
batlTC term, feature of asymptotic freedom, emphasized 

by lane’°, allows as to determine the large-momentum 


behavior of the complete kernel and consequently the large - 
“omentum behavior of Sj,(p^)> pmvlded S^(p^) does not vanish 
Identically, t^e approximation of sin^e-gluon exchange 

tn the kernel, Equation (X10) has two solutions (derived 
in Appendix which in the landau gauge are 


In'^CpW) 


33-2n^ (irregular) 

(3.11a) 

(tegular) 

(3.11b) 


regular solution corresponds to typical bound-state 
behavior, while the irregular solution corresponds to an 
elementary point-lifce ground-state pseudoscalar meson, 
Bi^dshy and Parrar^^ have shown that simple power counting 
in the vector-gluon theory with the regular solution 

implies the pion electromagnetic form-factor ^ 

asymptotically large (%(^^'" minimum 

number of fields n^ in the hadron). This is also the 
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experimentally observed behavior. For the irregular solution 
2 

one gets )^0('1)* Thus the usual experimentally observed 

bound -state physics of falling form factors, transverse 

momentum distributions, etc. can not be maintained with 

4-3 

the irregular solution. So, following Pagels , we conclude 

that the experiment favors the regular solution. In (3.11b), 

which is a constant that sets the scale for DSB, will 

47 

be the definition of the dynamical quark mass . Using 
the asymptotic form (p )~4m^/p , the amplitude 
below the strangeness threshold in e e"^ annihila.tion has been 
fitted with the conclusion that m^j ~ 244 MeV. According to 
a better estimate, using the logarithmic part of (3.11b) as 
well^®, m^ Si 315 MeV. 

As noted in. the introduction of this chapter, one 
can also proceed directly via the SD equation for the quark 
propagator and find out a XSB solution in a confining QCD. 

By confining QCD we mean one in which the effective gluon 
propagator could be written as in QBD but with the replace- 
ment (modulo gauge terms) 

2 j 

^\i.V ^2 "Sp y ^2 '2 

2 2 2 2 4. 

where g (k ) ‘^1/k as k -*• 0 arid C2 = ^ is the quark Casimir 

eigenvalue. This singular form, in a rough sense, corresponds 
to a linearly rising potential at large distances in configu- 
ration space. This form of the gluon propagator has been 



In particular, Baker 
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advocated by several authors ’ 

49 

et al, have found that in the infrared region in axial 
gauges singular part of the gluon propagator has the same 
spin structure as that of the free propagator and (the 
unrenormalized one in euclidean space) is given by 

jj(s) , ' Z(m) + MfL.) (3 12) 

where Z(m) is the renormalization constant for the gluon 
wave function at the renormalization point M and A is some 

35 

other constant. Ball and Zachariasen (BZ ) have used 
Equation (3.12) in the SB equation 

S“^(p) = S””'(p) + g^C2^(dk)Y^ (k) S(p-k)ra/ (p-k,p) 

(3.-13) 

satisfied by (unrenormalized) axial gauge quark propagator 

S(p) = B(p^) + G(p^) + (|(H + 1 ) 4 (3.14) 

2 

whereas the free propagator is Sq(p) = 2^/p • 5'or 
S(p-k)ry (p-k,p)S(p), BZ write an expression which is its 
longitudinal part chosen in a way that is consistent with 
the WT identity (which, in axial gauge, has the same form 
as in Abelian gauge theories); 
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, V _ (L) / N o/ \ r/l'+S” N +P )y y 

8(P' )r i; ^ (P' »P) S(p) = - [(^^^)Yy+ 2 1 

P' - P 


+ 


G' - G 

,2 2 

P' “P 


(fvl^+j^'yv')'^+ terms 


(3.15) 


Here p' = P“kj and F =F(p), F' =F(p')> etc. have been 
used. The equation thus obtained from Equation (3.13) can 
be broken into two equations which are equations for F and 
G separately if we choose the gauge p.n=0. We have taken 
this opportunity to work out the derivation of F and G (given 
in Appendix G), and found out that the correct expressions 
for the (renormalized) functions F and G are slightly different 
from the ones found by BZ; 


5'(p^)=- -V^^(1;1; - P^/PM^) (3.16a) 

PM 

G(p^) - i(2; 3/2; ) (3.16b) 

2 2 

where p= AOgZ (M)g^/4'Ji: . According to BZ, functions H and 

I may be dropped in case one is working in the gauge p.n=0. 

Henceforward, we shall assume this special choice of gauge 

for convenience, wherein we may set H = I = 0. To find 

out , instead of following BZ, we shall follow the 

1 fi 

normalization procedure introduced by Cornwall (remembering 
that our Green's functions are euclidean); 


C 


-1 

1 


S"'”'(0) 300 MeV 


(3.17) 
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If we take the asymptotic behavior of (3.14) using 
(3.16) and compare it with, the usual euclidean fora of 

fermion propagator for asymptotically large momenta, then 

^2 2 

we find that m(p )(rl/p which may be compared to (3.11b). 
Thus we see that the quark propagator defined by Equations 
(3.14) and (3.16) has the correct momentum dependence even 
for asymptotically large momenta. 


3.3 Dynamical Quarks Mass and Hadron Physics 

As stated in the beginning of this chapter, nonpertur- 

bative dynamical quark mass plays important role in low 

energy hadron physics. This has been shown in quite a few 

works. Thus there are calculations of pion decay constant^ 

47 

and electromagnetic fora factor of pion More recently 

it has been used for the resolution of the U(1) problem in 

o 53 

QOD and calculation of the n -►yY amplitude . The results ■ 
thus obtained are encouraging. We shall illustrate it by 
calculating the pion decay constant^ 

The axial-vector current for quarks is defined by 
Equation (3.3b) 

(x) = ^(:x)^^N^Y 5 (3.3b) 

with a trace over color indices understood. In the following 

we shall set renormalization constants = 1, and also 

2 a, 

A(p ) = 1 in Equation (3.4). \ 's are now isospin matrices. 
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Pion decay constant appearing in Equation (5.7) is 
giyen by 

J^5(0)I0> = -i (3.18) 

while quark -quark-pion vertex ,fp(p,p+l) Y 5 given by 

e^P^< 3 T®'(q) I T(f(x) f(0)X^)l0> = iS(p) (?(p,p+q)y ^ iS(p+q) 6 ^^ 

(3.19) 

Trace over color indices is understood. When p-integral of 

( 3 . 19 ) is substituted in (3.18), one readily deduces 

iq^^ 6 ^^ =/(dp) Tr[iS(p) '|^(p,p+q) Y 5 iS(p+q) Y 5 3 

(3.20) 

where the trace is over Dirac, flavor and color indices, and on 
the right the limit is to be taken as q '*0. Next, we expand 
r.h.s, of (3.20) retaining only linear, terms in q. With the 
approximation that the regular part of axial -vector vertex 
is Y.. Yc ^ /2, one has 

f„ — (P (k,k+q) I Y c = ""Yc 92^8 (p)‘"Y 2 ;Yc. (3.21) 

When Equations (3.21) and (3.8) are substituted in r.h.s. of 
Equation (3.20) in linear terms in q, it takes the form 

4=-12iY(dk) j- ^Sj)(k^)]. 

( 3 . 22 ) 



This last equation has been used in Eeferences [l6,4'7], and 

in a slightly modified form in Reference [ 53] to evaluate f ^ 

In Reference [47], S-p.(p with m^^ 244 keV has been used 

P 

with the result that f^ = 83 MeV,. close to the experimental 
value of 93 MeV, Similar results have been obtained in the 
other two references using slightly modified forms of 2^ 
to take into account its infrared behavior. 



Chapter 4 


ANOMAIOUS MAGNETIC MCMENTS OP IIGHT QUARKS AND 

dynai^iicaij simmbthy breaking 


4, 1 Introduction 

The magnetic moments of the proton and neutron are 
usually calculated in terms of confined constituent -quarks. 

A tetter knowledge of the magnetic moments of these baryons 
requires a better knowledge of the anomalous magnetic 
moments (a.m.m, 's) of their constituent quarks. As noted 
in Chapter 3, dynamically generated quark mass plays an 
important role in low energy hadron physics. In view of 
this fact, it will be interesting to study the effect of 
the introduction of dynamical quark mass on the a.ra.ra.'s 
of light quarks u and d, which we can assume, to a good 
approximation, to haTe vanishing current quark masses. The 
content of this chapter has been taken from Reference [54], 

In QED, the a.m.m. terra (in the lowest order) has 

the form 2 ^ 2 ^'" • In QCD where chiral symmetry 

has been broken dynamically, one naively expects the same 

form except for the change m -* and a -*■ But 

the dynamically generated mass depends upon the coupling 

constant in a nonanalytic manner such that it vanishes in 

10 11 

the limit of vanishing coupling ’ . That is, even for 

small couplings, we expect a.m.m, term in such a theory to be 


4 ^ 


pretty laxge. We shall illustrate this by an example. 

Consider a modified form of Nambu-Jona-Iasinio model 
given by the Lagrangian 

i (i^- lA(3^Ay -03, A,,)^+g[(i^ )^+( f iY^'SO^] 

(4.1) 

which is invariant under ordinary local phase transformation 

’?(x) -*• exp(i©(a:)) W(x), ■* 0(x)/e 

as well as chiral (global) transformation 

— — 4Yc0 

(x) ■+ e ^i'(x),, li'(x)-*' \?(x) e ♦ 


Gauge terms will be considered as a small perturbation 
over the original field theory. After the chiral symmetry 
has been broken dynamically, the fermion develops a mass 
S (p) = m given by^ (neglecting the effect of the gauge 
interaction) 

= 1 - ^ In (J^l + 1) (4. -2) 

S A ® 

2 

Pit 

where A is an ultraviolet cut-off and 0 < — ^ < 1 for 

, gA 

getting real A/m,. Next, we shall use the resultant fermion 
propagator 


s(p) = / Afp) 


(4.3) 
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to calculate the lowest order vertex diagram given in 
Figure 1, This gives the a.m,m. term; 


r ii(p,p+q) -Vji = 2igJ-((ik) S (k+(a ) 


iV 1*' 

2m 



SA 


r 2 ^^ 

L ^ 

gA. 



2 : 
A + m 


(4«4) 


2 2 

If 2n / gf\ is chosen very small then a.m.m, indeed becomes 
very large* 


In the literature various forms of 'running' 
(dynamical) quark masses have been derived and used for 
calculating low-energy hadronic quantities. We shall use 
them in the quark propagator in the following SB equation 
satisfied by the quark-electromagnetic-current vertex 
function 


rvi(p,P+a )=22 Y|_^ -/(dk) [S(k) (k,k+q)S (k+q)J K(p,q,k) . 

(4.5) 

QED will also be assumed to be operative but its effect on 
the mass splitting of u and d quarks, and its contributions 
to K will be ignored. 
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4.2 Solution3 Based on Asymptotic Properties of QCD 

In this section we shall first use the asymptotic 

solution (3.11b) (with Sp -♦ E and ignoring the 

2 

logai'ithmic factor. Also A(p ) in Equation (3.4) will be 

taken to be unity. This kind of solution has been used by 

47 

Bagels and Stokar to calculate the pion decay constant 
(as shown in the last part of Chapter 3) and the electro- 
magnetic form factor of the pion successfully, and also to 
calculate the quark electromagnetic self-energy, A consistent 
way to use this solution for our purpose would be to 
calculate the lowest order vertex diagram given in 
Figure 2. When this is done with the replacement E ■* m 
in the denominator, it turns out that the Pauli form factor 
F 2 ( 0 ) defined by 

rii(p+'i,p) = + gi-V i"' 

becomes complex. Evidently the asymptotic form of E , 

which we have used, raimicks a massless scalar propagator 

and is not suitable in the low energy region of integration. 

1 6 

As a remedy, following Cornwall , we shall parametrize 
S (p^ ) as 

S(p2) ^ . (4.7) 

•„2 . 2 


This gives the following expression: 
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Fig. 2 
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2 ^ 

g Cp m 


1 1 

^ dx J" dy [ — 

0 m^-q^yd-y) 


0 


2m (1-x) 


2 2 2 2 
m X + (A - m )(l-y)“q xy(l-y) 


^ (4.8) 


which, for q =0, reduces to 


g C 2 


^2(0) = — 21= 


+ 


^°2 


T{(/^ - a) 


m 


X 


[ In 


2 A A 

~T~7 T 

A A. “ni 


A^ 

T~^ + T 

m 


1 ] 


^ f (aW ) 


(4.9a) 


(4.9b) 


where m = 2 (p'^)j = " '~ 2^ ~ 2 '* provided that it can be 

l^=«i A -in 

solved for m for a given M and A . 

The first te^nn on r,h,s. of (4.Sa) is the standard 
result whereas the second term becomes complex for m >J\ . 

2 2 

The plot of f vs. A /m (Figure 3) shows that it 

2 2 

starts with a negative value and at A /m si 7.8, it 

2 2 

becomes zero. As A /m becomes infinitely large, it tends 

1 6 

to 1 / 2 , the standard result, Cornwall has given an estimate 

of M - 300 MeV and A 600 MeV, But with this set of 

values, a solution to the equation m,= —n — w does not 

A -m 

exist. As an example, the following set of values satisfy 
this equation: M 300 MeV, A ~ 900 MeV and m ^i360 MeVj 
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2 2 B S 

and for this set t{/\ /m ) -0,11. For the accepted value^^ 

2 

ttg (1 G-eV )?^?0.50 (and talcing maximum mass scale /\ above 
as the argument of in (4.9b)), one finds P2(0)~“l/45. 

4.3 Solutions Based on Infrared Properties of QCI) 

In this section we shall use a singular effective 

form of the gluon propagator. For the XSB quark propagator, 

we shall use either a parametric form given by Equation (4.7) 

or a form consistently derived using such a singular gluon 

propagator. Furthermore, following Acharya and Narayana- 
36 

Swaray , we shall assume for small momenta, where is 

expected to dominate over Kp that the leading (most singular ) 

terra in K has the ladder form with a singular effective gluon 

propagator, (k) k "This approximation to the 

dynamics may be regarded as an effective strong coupling 

"1 ^ 5 "t 

approximation". Works of Cornwall and Eichardson^ support 
this hypothesis. 

The form of gluon propagator introduced by Baker 

et al."^^ (Equation (3.12)) and the resultant XSB quark 

35 

propagator derived by Ball and Zachari.asen"'^ (BZ) (corrected 
form Equations (3.16)) is the one which we sh^l use for this 
calculation. Again we shall evaluate the one-gluon-ex change 
graph of Figure 2, but here the quark-photon vertex will be 
the complete one (all the calculations below are in euclidean 
motric); 
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,<:r. Nc 


/. /. 7., 

~1T!IT30 


ru(P»P+q.) = + go^2^ (dk)YxS(k)rii (k,k+q)S(k+qX£j WP“k)- 

( 4 . 10 ) 

Here we shall approximate the full quark-photon vertex by 

its longitudinal part. Since the quark -phot on and the quaxk- 

gluon vertices obey similar WT identities, following 

Reference [ 36 ] , identical longitudinal decomposition can be 

written for both the vertices (except for the difference that 

the quark-gluon vertex will have a color matrix), at least 

35 

for small momentum transfer. Thus, we shall write 


S(k)r^^^^ (k,k') S(k’) = - [^(P+F* )Y|i+ ^ (F’-F) x 


2 JjfY^J^' + (k^+k'^)Y^ 

^ ^ 

kf - k 

+ ( 0 »-G) -^—5 2—^ ] 

k* - k 


(4.11) 


for the quark-photon vertex as well for small q = k'-k. 

It is to be noted that if the mass-shell condition is not 
used (which is the case for the confined quarks we are 
dealing with), then the contribution to the Pauli form factor 
can come only from the last term in the square bracket of 
Equation (4.11). Thus we are interested in calculating the 
following integral 
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= / 


■2s -.,2 


(k-p) 


& te . ! ) - |(k } y^(y^yif + Y^ x 

k' - k 


[6 


7\a 


2 

(k-p)7^ncr+- (k-p)<^^ (k-p X (k-^) )cy n 

-f n J , 


(k-p) .n 


( (k-p ).n)' 


(4.12) 


Details of the angular integration have been worked out 
ip. Appendix H. Choosing all euclidean y antiher- 

mitian, the result of the integration is 


\ = q. y [Yj^jY^Jt^ 1(1 + ^)(2 - y * ) / 


2 , (p+q) ,, 2 


dk‘ 


(p+q)' 


2 2 


(k +p +p.q) 


qV4 


,4 ,, 2 2 . 1/2 

(q^ - 4k q ) ^■ 


- J dk G(k ) [ " ■ 2— - g 


^ ^ /s ^ ^ ^ / 4 2 2sl/2 2 

0 q (k +p +p.q) (q -4k q ) ^ (k +p +p„q) 


q terms + n.. terms 

p p 


(4..13) 


x^rhere the constant vector n has been chosen to be perpen- 
■^icular to the hyperplane defined "by p,q and (icr^^yq^,). We 
have not been able to locate any formula in the literature 
which can be used to perform various integrations, which 
appear in (4.13), analytically, if the transcendental form 
given by Equation (3.16b) is used for G(k ). However, the 
integration can be done in certain limiting eases: 
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2 2 2 

(i) When p ,q << .pM so that the series expansion for 
function is valid and it can be approx’imated by the 
first few terms in the series; then 


.^)(2 

q (p+q) pM 


+ 4 )(— 1 


1 


P 


) 4(q +2p.a ) 


p.(2p+q) (p4-q).(2p+q) 5P M' 


(4 ^ 8 ij"+e^<i 1) In l atq)- . (2p-Hi) ] - _1_ 

3pM 5p 1-1 p.(2p+q) ■' q^r 


— [2 - (r+ 1 
3 pM 


in 

(r+1)^/ -1 


i -4-K - 2r (r+1 In . ] } 

5 o2,,4 h (3 ,+t)V^_-1 


p 


+ q terms + n terms 


(4.14a) 


where r = ^ . Now on taking the p 0 limit: 


-qy [y^, Y I C 2 Q 


(1 — ^ + q terms + n,, terms 

r- r;'. II. LI 


3P>? ' 5pM" ■ '1^ 


471 ? 


^ - l 2 >tYa; , Yn] C. + q terms + n, terms. 

^ ' b* 1 -tQvA U U. 


q -♦ 0 
This gives us 


3 Pm 


(4.14b) 

(4.14c) 


1 4 


r p.(0,q) = 3 ,‘*■'^11 


(4.15) 



wiiere dots indicate other possible lorenta structures. 


(ii) When q - 0 (but with finite p): 

P 


2 

&(P^)>+ q^_^te3nns 


n terms 


(4.16) 


This gives, 


rp(p,p+<l) = - G.^Ci/2-S(2;3/2;-^-^) J 

+ Y,, terms + (4.1 7a) 

If we substitute numerical values fcsr the quantities 
appearing in liquation (4,l7a), namely,^ (M^ )M^ « x 
(0.16 GeV)^ and p^Ji(0.5 GeV)^ then we get 


Tn (p,p+q) b; - 


i ^\iv 1 y 1 


2x300 MeV 


+ Y,, terras + . . . . 


9 ■ ‘P 


(4.17b) 


Similarly, 

ru(P+9,p) 

^ 2x300 MeV 



teitns . . 


(4.17c) 


In minkowskian metric, this can be written as 


TpCp+GiP)/ 


i cy, 




'\i.v 


1 mink 


2x300 MeV 


X terms + 


(4.17d) 
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where we have followed the common practice in 

^ 2 

assuming the same form for nijjyjjCP ) in the spacelike 



2 

region (p < 0) as in the timelike region. 

We have taken the same form (Equation (4.11)) for 
the longitudinal part of the both quark-gluon and quark- 
photon vertices. Assuming that their transverse parts 
are not much different from each other, we can write 
the counterpart of Equation (4.17d) for the quark-gluon 
case as 

_ n 

ru(P+<l>P)-" — X § +Y terms +.. .. (4.17e) 

^ 2x300 MeY ^ ^ 

where C is a number, likely to be 0(1 ), which takes 

into account the non-Abelian character of the gluon. 

There will be a similar change in Equation (4.15). 


(iii) When q^>> p^*^ p]/ : Since momentum transfer is no 
longer small, strictly speaking relation (4.11) can 
be applied only for the quark -gluon vert^.. We are 
assuming that when Equations (4.11), (4.12) and (4.13) 
aro used in Equation (4.10), the result, so obtained 
can be extrapolated to large q even for the quark-photon 
vertex . Thus we shall approximate the quark-photon 
vertex by 


r^(p,p+q)2: 


SL 




4p q 


+ Y„ terms (4.18a) 

H' 


f p^(P»P+<l) ^ 


d..„q^ Li^+v terms +. . . , 


.22 ’P 

4p q 


(4-l8b) 



r li (p+q,p) 


(4.18c) 


+Yn terms +. . . . 


r’uCp+qjP) 


mink 


4p^q2 ^ 

,2,, ,4 

.. |: 

lo^^ioq 


ivrc, 

+Yjj, terms+. 


a. P 


4p q 


. (4.18d) 


Similarly, we assume without explicitly calculating 
the triple gluon vertex contribution, that the quark- 
gluon vertex has the same form [as Equations (4.18a)- 
(4.18d)] apart from a constant factor, O’, 0(1). 


If we use expression (4.7) for S (p'^ ) (and the 


mass-shell condition) along with the singular part 

,2 ,2 

of gluon propagator D j,(k)-^ g,.i; ^ S’ » 

^ ^ ^ (k -1^ +ie) 

which has been used by Cornwall’® for calculating f^ etc, 

then the Pauli form factor becomes linearly divergent 

as p. -*0. Cornwall has argued that p should be 

kept finite and equal to a typical hadronic mass scale. 

On the other hand if, following Cornwall, quarks are 

assumed to be confined particles having no mass shell, 

then the Pauli term does not arise. 
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Acharya and Narayana-Swamy have also examined 
the feasibility of dynamically broken chiral symmetry 
in QGD with zero -bare-mass quarks interacting via 
single gluon exchange when the gluon propagator has 
the infrared behavior; D^y(k)«'» k They have found 
that the chiral symmetry can be, although not necessarily, 
realized in the Nambu-Goldstone mode with 



(4.19) 


—1 2 

S (p) = -m(p ) = constant . 

As can be easily checked, this result is too strong to give 
the Pauli form factor. 



In this section we shall explore some of the experi- 
mental consequences of the presence •f. Pauli terms given 
by Equations (4.17) and (4.18), 


(a) Baryon magnetic moment; If we make the usual assum- 
ption that the baryon moments arise solely from 

the constituent -quark moments, then following 
57 

Barik and Das , and references given therein, 
we can obtain expressions for the magnetic 
moments of the proton and neutron in terms of the 
magnetic moments of the corresponding constituent 
quarks. This has been done by writing the baryon 
magnetic moment as 

113 = S<Bti cj (Bt>. (4.20) 

where I Bt> stands for the state vectors of the 
baryon in question and in the present case, it 
represents the regular SU( 6 ) state vectors. The 
well known relations between the baryon magnetic 
moments and the etwrresponding constituent-quark 



moments are; 


= 5 - ^‘d). •‘n = T - ^u)- (+- 21 ) 

If we take a simple expression for the quark 
(Dirac) magnetic moment in units of the nuclear 
magneton as 

D ^ 

®q (4.22) 


where Hp is the proton mass, ra^, the effective 
mass of the quark defined by Equation ( 3 , 17 ), 
and e^, the electric charge of the quark in the 
unit of the proton charge. Thus from Equations 
(4,22) and (4.l7d) we have (in nuclear magnetons) 





1 / 9 )- 


( 4 . 23 a)' 


1J.^ 5 :^ 1 .86, ri - 0.93 . 


( 4 . 23 b) 


This result, when substituted in Equation (4.21), 
gives 

Ppii 2.79 , -1.86 ( 4 , 24 ) 

which can be compared with the experimental 
numbers 2.79 and respectively. The agreement 

may be fortuitous, particularly in view of the fact 
that the value chosen for m was a bit uncertain, 



but what is remarkable is the fact that the 
dynamical mass, given by Equation (3.17)(which 
was so chosen in a different context, namely, 

*1 ^ 

in the calculation of the pion decay constant ), 
together with a.m.m. given by ( 4 . 1 7 d) (where again 
the same dynamical mass was used) can yield a 
number which is so close to the experimental one, 

(b) Spin-dependent potential energy between quarks: 

Here we shall calculate the spin-dependent part 
of the potential energy between a heavy quark 
and a light antiquark. The light antiquark 
may be regarded as moving within the color 
field provided by the heavy quark. Our argument 
below will not be as airtight as we might 
wish. For the sake of simplicity, we shall 
depart from our earlier convention and assume 
that the quarks may be treated as on-mass-shell 

■jf 

within the framework of a confining potential. 

^ArbuzHov^^ has shown that quarks may have a mass-shell even 
though the gluon propagator of lA~type is used. Although 
m* calculation of the Pauli term does not take this possib i lity 
°nto account, we shall do it here and in the following due to 
our igmorance of methodology in treating such cases otherwise. 



w V-/ 


^ simple calculation of the spin-dependent 
part of the potential with one -gluon exchange 
in the no n-r el at i vis tic limit shows that 


S 


I fr -0 


2r <yi»PiXpJ 

& l± --' - - j— (1+2K^) 

4m ^ 




2m^m2 


4rar 


iop.p.xp' 

■ ' (I+K^) 

2m^in2 

4m^m2 


(1+K^ )(1-«2) ] (4.25a) 


whei*® and ^2 Pauli matrices which act 
on particle 1 and 2 spin wave functions respectively 
and Kp are the (constant) cbromomagnetic Pauli 
form factors and q is the momentum transfer; 
q = p^’-p^. Following the conventional wisdom 
that 1/Q.^'’‘type of potential corresponds to a 

2 j_jjearly rising potential in the configuration 
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space 


, we can Fourier -^transform Equation 


(4.25a) as 
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b- and c- quarks have large current-quark masses 
and hence and are, presumably, small; 
hence (4.26c) is essentially the result obtained 
by Eichten and Peinberg^^. This result agrees 
well with the experimental numbers. However, 
Equation (4.25b) has fewer free parameters than 
Equation (6,2) of Eichten and Peinberg. 

(c) Reaction cross-section of pp hadrons: Here, 

we shall calculate essentially the reaction cross- 
section of qq -*■ qq through one-gluon exchange 
(Figures 4a and 4b), For quark-gluon vertices 
we shall use Equations (4,18), along with the 
constant factor O’, and see the effect of the 
Pauli terra over and above the bare vertex. Since 
here wo are looking at the processes involving 
large momentum transfer, we shall use the 
conventional form of the gluon propagator 

• External quark and atttiquark lines 
are assumed to be on mass-shell. If 6 be the 
scattering angle and (p+k) = (2E) , then the 

differential cross-section for this process in the 
center of mass frame (retaining only linear terms . 
in Pauli form factor) is given by 
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4 2 /q 1 + Cos^e/2 2 Cos^0/2 


K £4Ji [ 

1287t^E^ Sin^ &/2 


Sin^e/2 


+ 1 (1+cos^e) ] + ■ ^— 

128^E 




mS 


l4Cos^©/2 QoB^Qr/2 003“^ 0/2 ^ ^ 

2Sin^e/2 2Sin'^0/2 2Sin^e/2 ^ - 

^4 2 /^ G'P^(M^)m'^C-, ^1+Cose 

128?? m? 2Sin^0/2 

6+Si?0+2Cose S±n9 + 2 Cose 

+ — 7' — + 5 J 

l6Sin^0/2 8 Sin 0/2 

(4.27) 


where the first two terms are the contributions 
ooraing from the bare vertex part and the third 
term is the contribution from the Pauli term. 

It is clear that for high energies, the contri- 
bution of the third term can be significant only 
for very small values of 0. To get an estimate 
of the angle at which the contribution of the 
third term is of the same order or one order of 
magnitude smaller than that of the first term, 
we have 


C»P^ Ci 


1 


mE 


2sin'^0/2 


1 or 10 


(4.28a) 



Sven for m ti 10 MeV (a number which is usually 
assigned at high energies) and GeV, this 

gives 

X 10"^ or 10“^. (4.28b) 

Even if we use the full- Dirac form factor (instead 
of just unity), the crude estimate (4.28b) is 
unlikely to change. Thus it will be difficult 
to detect the effect of the Pauli term this way. 

If, instead, we had looked at the differential 
cross-section for e qq, then in the equation 

corresponding to Equation (4.28a), there wauld 
not have been any (Sin0/2) factor on l.h.s, 
and the effect of the Pauli term would have been 
insignificant at all angles. 

4.5 Summary and Conclusion 

We have seen that anomalous magnetic moments (of the 
fermion s)can be significant, although not appreciably large 
as speculated in the beginning of this chapter (at least for 
the realistic cases), in theories where mass generation 
occurs through dynamical symmetry breaking. To calculate 
a.m’.ra.'s of li^t quarks we essentially used one-gluon- 
exchanga approximation ard employed solutions of both 
kinds; solutions based on asymptotic and infrared properties 
of QCD derived and used by various authors. We found that 




all solutions can not be used for this calculation. Solution 
-for the gluon propagator derived by Baker et.al."^^ and the 
resultant chiral symmetry breaking and confining solution 
for the quark propagator derived by Ball and Zachariasen^^ 
form a good combination for this calculation. The result 
thus arrived at vas used to check the experimental consequences 
of the presence of a.m.m. (both electromagnetic as well as 
chromomagnetic ) terra. It gave very good results for proton 
and neutron magnetic moments when regular SU(6) state vectors 
for baryons were used. In fact the use of Equation ( 4 . 21 ) 
alone gives the correct ratio of magnetic moments of proton 
and neutron, using dynamical mass (here the effective quark 
mass) given by Equation (3.17) and the resulting a.m.m, given 
by Equation (4.l7d), one obtains correct values for the 
individual baryons. In the case of the spin-dependent potential 
calculation, we obtained a result which agreed with an 
earlier calculation (which had a good agreement with experi- 
mental results) and, moreover, had fewer parameters. In 
differential cross-section calculation of qq •+ qq, we found 

that the effect of the presence of the Pauli terra was 

-2 

insignificant (at high energy) at any angle 10 , while in 

case of e^e'^ -» qq the effect of the Pauli term was insignificant 
everywhere. 

Although we have worked in the axial gauge (with a 
specially chosen direction of the axial vector n ) , our 
expression for a.m.m, is apparently free of n . 



Chapter 5 


DP- Al\iD DOWN- QUABK MASS DII'PERSNCS AND 
DYNiUyilCAL SYMMSTRY BREAKING 

Encouraged by the applications of the dynamical 
quark mass as illustrated in Chapter 3 and the results 
derived in Chapter 4, we shall use the dynamical quark mass 
to calculate u-d quark mass difference. This is the modern 
version of the old problem of p-n mass difference. If 
the origin of the quark mass lies in the dynamical breaking 
of the chiral symmetry of the QGD Bagrangian, then this 
may provide a mechanism through which QCD plays a role in 
determination of u-d quark mass difference (taken to be 
of electromagnetic origin) through its feed-back on electro- 
magnetic self-energy. With a new interaction coming into 
picture, one may hope to get the correct sign of u-d quark 
mass difference. 

2 

let us substitute Equation (3.4) with A(p ) = 1, 
in the gap Equation (3.10). Then one gets (dropping the 
subscript D in Sj,) 5 

Y 

v.2(p^) =-f(ak)C g 5 jK(p,i,,o) . (5,1) 

^ . k - S (k ) - 

With K as the QCD kernel for qq scattering, for large 
momenta in the approximation of single-gluon-exchange, 
the solution of Equation (5.1) was found to be (3.11b). 
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To proceed further, we follow Ref erence [61 j . 

Consider the self— energy of a quark obeying the 
integral Equation ( 5 , 1 ) where K includes both QCD and 
QED interactions. We assume the QED interaction to be a 
perturbation to the QGD interaction, i.e,, that K = K 

7 7 Q 7 

where is the QCD kernel in the lowest order and 
6K = , Here is explicitly the QSD kernel 

for qq scattering. is the effect on the kernel due 

to 6E , the change in the mass operator. Again we shall 
approximate 6 k by the lowest-order term in The 

lowest-order term in 6X^2 has been neglected because we 
have already neglected the two-gluon-exchange contribution 
to Kq. Writing E= 2 ^ + 62 we get the following integral 
equation for 62 : 


k^ + 2^ (k^) 


62 (p ) =/(dk)[ Y56E (k^) p 2 r-. T ' Y Z ]Ko(Ptk>0) 

Lk S- (k 


Yc 2 (k^) 

/(dk)[ I * 

k - S^(k ) 


( 5 . 2 ) 


The second term which is the driving term in the above 
equation is exclusively 62 After performing a Wick 
rotation and doing the angular integration, Equation ( 5 . 2 ) 
can be written (in the landau gauge) as 
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(p^) 


3 Cg 
16 # 


P 

/ 

0 


ggCk^) [k^ - 2 2 

« — „ ^5-^— — 62 )k^ dk 

P [k+So'’^)J 


+ 


3 Cg 

T?? 

•Z 2 

3 e . 


00 

s 


p 

/ 


2 2 [k^ - 2^ (k^)] 2 9 

ggCk ) s — 2 . (k'^)dk 

[k^ + ni (k^ ) r ^ 

2 2 2 2 2 

Sj (k ) k dk 3 

7'[k2+ S^(k^) ] le”? "" 


00 

/ 

_2 


(k^) dk^ 

? + S^ (k^) 


(5.3) 


2 2 2 2 

where we have let gg((k*-p) ) -♦ gg(k ) and i=(u,d). For 
2 

large p , IBquation (5.3) can he converted into a differen- 

2 2 2 

tial equation (upto terms of order T q(p )/p which we 
will ignore) of the form: 


[p^esi (p^)]" 


62 (P^) 

• Y 

p In p /A 


3a:. 4 m:!, 2 y 

(p ) A ■ 


4 % 


12 


33-2Nt; 


(5.4) 


2 

Here the differentiation is with respect to p and we 

2 2 

have used relation (3.11h) replacing ti by j\^ and 
2 2 k^ 

g (k )/4 TE =T?Y / In has been used. A is the usual 

® A 

free mass scale parameter of QGD, Equation (5.4) can be 
solved analytically. The solution satisfying the appropriate 
boundary condition 
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lim [ 6 S . (p^ ) ] ' ^ 


0 


2 

P - 00 


is given by 


6\(P-) 2 

P -► 00 


4 m- 


?r(UY)r{2n) J [fTsIff^ 


4 % 


n=0 


(In pVa^/+'' i- (n+2) 2 n+1 

][ In Eb- ] 


r(Y+ n + 2)f(Y+ n + 3) 


(5.5a) 


3^ . 4 

4# — 2^rO+Y)r(2+V)Cr(1-r)ii'i(1-Y;2ilnpVA® 

!P 


(In pW f 


r(Y+ 2)r(Y+ 3) 

.3 


5, (1,2iY+2,Y+3!ln pVa^)^ 


2^2 


2/2 

In p /A , 


(5.5b) 


3 a, '4 mi. p o 

nf r(1+Y)r(2+Y) ^(pVa^) (5.5c) 

P 


where are generalized hypergeometric functions. 

We observe that (2^ "’^d^ will have the correct 

2 2 

sign only if ^(p /a ) is negative. , Numerical computations 

2 2 

show that this does not happen in the range 1<lnp/^ <18 

2 2 

for Np = 6 and in the range 1 < In p //\ <16 for Wj, = 7 
(see ELgure 5). If A f® taJten to be 150 MeY, then this 
corresponds to (0.4 GeV)^< p^< (1410 Gelf' for Nj, = 6 and 
(0.4 GeY)^< p^< (520 GeV)^ for % = 7. 



(zV/jd)u7 
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Thus we see that for reasonable values of momenta 
we do not get the correct sign of ( E ^ 

in the case of Reference [47] > 62 . (p^) is completely finite 

2 

and tends to zero as p tends to infinity. 

At this stage, one may conclude that the persistence 
of the wrong, sign may be , as enumerated in Reference [ 47 ] , 
due to one or more of the following reasons; (i) There is 
a bare up- and down- quark mass difference, (ii) The result 
depends on the nonasymptotic part of 2|^(p ). (iii) Correct- 
ions higher than the lowest order (at least in the strong 
coupling) are important, (iv) Instanton contributions to 
the kernel are large. 

62 

However, recently Chang and li have also examined 
this problem and have found the correct sign for 
within the framework of with zero current- 

quark -mass and without invoking instantons. Their work 

r 'T 

is based on a new approach developed by Chang and Chang to 

calculate the dynamically generated quark mass in QCR 

g 

based on the Nambu-Jona-Iasinio mechanism, Chang and 
Chang have done a renormalization group (RG) analysis of 
the Nambu-Jona-Lasinio gap equation and found that it is 
indeed an RG invariant for QCD. With the renormalized 
two-point Green's function for the quarks given by 

,2 2x 

(p < < m ) 

p( 2 ) ^ 22”' {i - m) 


(5.6) 



thsy find that the mass of the quark so generated is^ to 
two loop RG accuracy, 


m = e1/6 


(5.7) 


(2 ) 

where is the free mass scale parameter of QCD evaluated 

up to two loop. This mass, m, is a EG invariant. Chang 
and la have shown the gauge independence of the result (5.7) 
too. To generate the u-d quark mass difference Chang and 
Id. include the effects of QED as we have done above. To 
one loop RG accuracy, the degeneracy between u and d quarks 
is not lifted. To two-loop RG accuracy they find 


m, 


u 


m^ ^ -ffl 


288 71 


'13 


(5.8) 


where d.j^ represents the genuine two loop constant of the 
self-energy graphs involving both gluon and photon exchanges 
and is approximately 64. This gives the right sign and also 
the right order of magnitude for 

The analysis of Chang and Iii has the advantage over 
ours in the sense that their result is RG and gauge invariant 
and is meant for the low energy region. Calculation up to 
two loop EG accuracy gives an additional advantage. 

In conclusion, we see that if QCD feedback on 
clectmmagnetic self-energy is properly taken into account, 
it may correctly account for the u-d quark mass difference. 



Chapt er 6 


C0NC1IIDIU& RBl^iiffiKS 


In the present thesis, we have studied two types of 
problems: (i) Dynamical symmetry breaking (DSB) in the 
Abelian gauge model of Jackiw and Johnson”*^ with the help 
of the gauge technique, and (ii) Applications of dynamical 
chiral symmetry breaking (XSB) in QGD, 

In the first case, the quantity of interest to be 
calculated was the axial -vector mass to the fermion mass 
ratio (M'/m). To apply the gauge technique to this problem 
one needs to know the fermion spectral function (P ). In 
view of the complexity of the problem of finding P , we 
chose to work in a sector where p/m was very small. If we 
assumo that the result for the mass ratio obtained by earlier 
authors is at least qualitatively true, then this implies 
working in a sector where the axial -vector to vector coupling 
constant ratio (g'/g) is very small. In such a case ;ve 
can hope to expand p in terms of these small parameters 
(unless P depends upon these ratios in a nonanalytic way). 
This seems sensible because, for example, we expect P to 
reduce to the QBD P smoothly in the limit g'/g - 0. Since 
the gauge technique is known to work well in the infrared 
region, precisely where DSB- is supposed to occur, we expect 
it to give reliable results even in the aforementioned sector. 
To our surprise, we find that this way the result for p/ra 



comes out to te extremely small. This does not conform to 
the results obtained by earlier authors in the above stated 
sector. Firstly?-, using the gauge technique, we do not agree 
with earlier authors that one can use axial -meson propagator 
with mass p, = 0 to obtain a non-zero fermion mass and then 
subsequently to obtain p / 0 (see Chapter 2). This seems 
natural, since unless one takes the perturb8.tion expansion 
around the new vacuum in which the particles have already 
become massive, conventional perturbation theory will not 
reveal the structure of the massive theory, if it is a 
massive one at all. Secondly, we evaluate p /m with a 
terra in P which contains information about p/ra being nonzero 
and thus p/m is determined in a self-consistent way. Our 
method cannot determine the individual masses (because the 
mass-scale in the theory is set essentially by the renor- 
malization process, which we are not using explicitly). In 
particular, wc can not say whether the mass ra is non— zero-. 

In the second case, we have applied (dynamical) XSB 
solutions for the quark propagator to two problems: 

(a) Anomalous magnetic moments (a.m, ra.'s) of light quarks 
and (b) Up— and down— quark mass difference. 

We have seen that in theories in which chiral 
symmetry breaks dynamically, a.m. m. 's of quarks can be 
significant, although not appreciably large as one might 
speculate naively (see Chapter 4), Among the various XSB 
solutions for the quark propagator that we have examined. 
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35 

the solution obtained by using the singular form of gluon 
49 

propagator and the longitudinal part of quark—gluon vertex 

is the most interesting one. To be consistent, we too have 

used the singula.r form of gluon propagator in the calculation. 

Since the quark propagator which we use is supposed to describe 

confined quarks (because it has no "mass shell" singularity), we 

have not used any mass-shell condition which would have 

otherwise given an additional contribution to a.m.m, term. 

One interesting aspect of our result for the a.m.m, term 

is that it has desirable momentum dependence in the infrared 

as well as in the ultraviolet region, i.e., it vanishes 

linearly with q (momentum transfer) as q - 0 and inversely 

with q as q -* 00 , This result (for low energy) along with 

the effective quark mass was used to calculate proton and 

neutron magnetic moments using regular SU(6) state vectors 

for baryons. The result thus obtained was in good agreement 

with experiment. We also calculated spin-dependent potential 

between a heavy quark and a light antiquark. When this was 

applied to calculate the spectrum of ground state mesons, 

we found that the anomalous chromomagnetic moment had no 

magor role. We further found that the contribution of the 

anoraalous chromomagnetic moment in differential cross section 

for qq -* qq was insignificant (at high energy) at any 
-2 

angle © > 10 . 
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Our result of the calculation of up- and down-quark 
mass difference comes out with the wrong sign for any 
reasonable value of momenta, although it is finite and tends 
to zero as the (momentum) oo , (see Chapter 5). We are 
left with one or more of the following reasons for the persis- 
tence of the wrong sign; (i) There is a bare up- and down- 
quark mass difference, (ii) The result depends on the 
nonasymptotic part of the quark self-energy function' . 

(iii) Corrections higher than the lowest order (at least 

in the strong coupling) are important* (iv) Instanton 

62 

contributions to the kernel are large, Chang and li 
have obtained the correct sign (and also the correct order 
of magnitude) of m^ - m^ from their calculation of masses 
at low energy up to two“ioop renormalization group accuracy. 
This work narrows down the possible ways to improve our 
result. It indicates that it may be possible to obtain 
“u ^ ®d generating the quark masses entirely 

by color and electromagnetic forces. We feel that a possible 
way to improve our result may be to analyze the XSB 
solution of the Schwinger -Dyson equation for the quark 
propagator in the infrared region using both gluon and 
photon exchanges, and using longitudinal vertices for 
the complete ones. The result may be further improved by 
requiring it to be gauge Independent and renormalization 
group invariant. 



78 


We have not shown the gauge independence of any of 
our results in the thesis. This is a difficult problem and 
only recently some progress seems to have been made in this 
direction. Our attempts to use the idea of DSB to extract 
physically observable results are encouraging to consider 
further work on it worthwhile. 
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Appendix A 


INTJilGRAL IQUATIOM FOR FFRffiON SPFCTRAX’ FUNCTION 


The more familiar form of _ spectral function given by 
00 + mpg (p^Pdp^ 

s(p) =/2 , 

m p -p + le 

with = Z2 / m P 2 (p^)dp^, 1 = Z2/ P^([.i^)d;p^ is connected 

22 

With the form (2.3) by 

P(W) = e(W)(WP^ (IV^) -f mP2(W^)). 


Now using Equations (2.4a), (2.4b) and the spectral 
decomposition for Z 2 in the SD equation for the fermion 
propagator (Equation (2. 12)), one gets 

where 

p ^ \Ji2^ 

\(p,w) = ig r (dk)Y^ 'fy K'i (k) (AH-S) 


2b(p,W) = ig' / 


(dk)Y^tY5 




-VS 




(A-li) 

Pole term in axial vertex does not contribute to 2^ (p>W) 
because gauge boson propagators are . in the landau gauge. : 

t ' ' . '' ■' ■' ' V 7'/ 

Using projection opei'a>tor#f ^ SW , which converts 

iPr.'vjV 4 . f 



80 


equation A(p ) + B(p^) = 0 into /p^ A(p^) + B(p^) =0, 
TSquation (A-1 ) can be written as 


/ dW'_P (V/M 


WWTlkw) -Sg(W,W')j= 0 . (A-2) 


On taking imaginary part of Equation (A*-2), it becomes 


7[e(W) (•' (W)r W - (W,W) - = x 


If we define 


m = 2^ (m,m) + Sg(m,m) , 


Im[s_^(W,W' )+S 2 (W,W' ) ], 

(A-3) 

(A-4) 


then we can write 

P(w) [w - (W,W) - Ej (W,W) ] . 

= P(W) L W-ni+S^(ra,m)-Sj^(W,W)+ EjC®*®) -*Sg(W,W)J. 

p(w) - S^(W,W)+ -2^ (W,W)] is 0(g'^,g^,g^g^) , 

i.e., of the same order as the transverse corrections which 
have not been taken into account in any case. Thus, upto 
these corrections, Equation (A-3) can be written as 


■e(W)P(W) fw-m] = Im[E^(W,W') +S 3 (W,W')]. 


(A-5) 
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Appendix B 

liVLlGINARY PART OP P:MION-S1 IiP-P,UERGY PUI^CTION 


Pepmion self^energy due to its interaction with 
B-meson in the landau gauge, to, the lowest order is 


Vp) ^ 


k^kj; 


V 




(2^)4 5 ff-X-m+ie ^ 


3m+j|( - 31^ + 


. ^ |2 r d^k 

— "f'Xg * i 

(PTT)^ [(p-k)^-m'^+iej(k'^-ii^+ie) 


k +i£ 


= p*(p^)+^( p(p^) 


(B-1) 


where 


)?-=ig,2/_lk ^ 

(27x)^ [(P“k) “tci +ie](k -fi^+ie) 


and 


p2 P(p^) 


A P^-3p.k + 2XPikl 

ig' 2 " 2 

(2'jc)^ [(P”k) -m +ie](k -p. 4-is) 


(B-1 a) 


(B-1h) 


Pirst consider the integral in (B-la) 






(270 (pQ-kQ-dffie)(pQ-kQ+A-ie)(kQ-K+ie )(kQ+i<-ie ) 

A = {(?-k)^ + h = (£^ (B-^a) 


We first perform the k^ -integration by closing the contour in 
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the upper-half-plane; 

^ / d\ r 2® 

(2^0^ 2A(p^-A+H)(Po'-A-K+ie) 


(p^-A+K )(Pq+A+k)(-2^) 

Equating imaginary parts on both sides ve get 

,2 „ 6(p-A-}() 

Im F’ = ^ / dQ diSl (B-2c) 


Ikl =p, 

(B-2d ) 

where 9 is the angle between p and and p is the solution 
of the equation -a-h = 0 tA I kl : 

2 ^ 4 2*2 2^^ 2_ , 4 2,1/2 

ft ^ iptCose + no +4|i PqP Cos e-4 p u ) , 

p “ p ir — 2 V 

2(Pq - P Cos e) 

^2 2 2 2 

X = P + - m . (B-2ej 

After a straight forward angular integration, Equation (B-2d) 






'M 

4 


/dQ 




2 Ik( Ikl - Ipf Oos0)+llcj Al 


becomes 



Im F' 





(B-2f) 


e-function on r.h.s. of (B-2f) insured that the expression is 
real. Integrations in (B-lb) can similarly be done; 


Im ( p' ) = Im ig'7 




(211)^ [ (p-k) -m^+iiE](k^-'[i^+ie) 


1671 ''\i 


(P,m) , 


(B-3) 


Ira (p'Fg) = Im ig' / 


. .2. -3 P.k 

(271)^ [ (p-k) -^n^+ie](k -p +ie) 


ImCp^F^) 


32 up ^ 


(B-4) 


Im i..2/ A^k- 2( , p,k )^. 


( 271 ;)'^ [(P'“k) ‘-m -fiG:^(k -P' +iE)(k +ie) 


[(p^+p^-m^)^^(p,m)-(p^-m^)^0Q(p,m)] , 


2 2x2, 


32 Tip 


(B-'5) 


where jZ^^Cp^m) = ;2^p,m)|^^Q. Thus, 


Im F == Im (F^ t Fg + ^ 3 ) 


(B-6) 


32 Tip 



Subatitutlng Bquations (B-2f) and (B-6) in Bquation (B-1), we 

get 

52 p __ (J’ 

+ ^ ^ 0^) J, (B~7) 

P 

where 

0 ,,, (P,W) = { (p^-(W+p)^ {p 2 _(y_ 1^)2 j 1/2 

= (P^-W^) e(p^-W^) . (B-8) 


Taking the limit 


lim \f {0^-- 0^) 


•{(P^+W^)/(P^-W^) } 0(p^-W^), 


ji-^O 
Thus, 

Im E. 


'B 


Uf^) ( 


|i”*0 


' ^ ■ -^1^ - (p^-W^) 0(pW). 

167ip 


(B-9) 


Now, from the definition 

V ' \ 42 , d\ ^ 1 

Sjip) = -1« 


^ 1 - U H 3 ; , 2 


k]^ 


le 


k +ie 

”2 

k + ie 


■) 


(B-10) 


it is cl oar that 

2 

Im SA(|f,W) = + (p^-w^) e(pW). 

16 7i:p^ 


(B~11) 
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Appendix C 

FSRMION -SPSCTRAL PUNCTIOIM lU QSD 

Integral equation satisfied by the fermion spectral 
function in QSD is Equation (2.16) with =0: 

2 

0 €(q) ((«>- 1 )r°’ («) = - 6 /do'r°^ (cjJ' ) o ' (o+w ) . (0-1) 

let us substitute 

e (o) ) r°’ (w) = w r^ (ci)^) + r2(w^) 

in Equation (C-1). On equating even and odd parts in w , 

2 

we get two equations (o) = z) 

z 

r2(z) -r^(z) = -8 / r2(z’)/z' dz', (C-2a) 

z 

z (z) - r2(z) = -6 <^ r^ (z* ) dz’ . (C-2b) 

This pair of equations can be reconverted into hypergeometric 
equations 

2 

[ z(1-z) —^ + {l-z(3.+2 6) } ~ - (1+6) ]r.(z)=0, (G-^3a) 

dz 

2 

ljz(1 -z) — ^ -2(1 + 6) z ^ (1 + ^) J 1*2 (z) = 0. (C-3b) 

d2i^ 

The appropriate solutions, satisfying the boundary conditions 
embodied in (C-2a) and (0-2 b) are 
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(z) 


ro(z) 


ilrl) 


- 1-2 


r “ 26 ) 2 - 2 ^ ?{-6,-6, - 26 , 1 - 2 ), (0-4a) 




(G-4b) 


Normalization in (C-4a) and (G-4b) is such that 


r° * ( w) -♦ 6 (CO- 1 ) as g ■» 0 . 


Thus we get expression (2.17). 

Now we shall show that the G-oldstone mode can not be 
realized with a fermion spectral function of the form r^’(co) , 
Let us substitute expression (2.21) in Equation (2,11b); 




// 




r 


a 2 


(C-5) 


where lira a 1 will be taken after the integration. We 
shall use the standard integrals 


•/“ P(a,b}c;-x)x dx 
0 


r ( a+3 )r (b+s )r' ( c )r (-si 

f (a)r(b)r(c+s) 


and 


/ x®** (x+y)"^ P(a,b;c;'-x)dx 
0 


rra-c4d)nt>-c+d)r(c) 

r(a+b-c+d)r(d) 


X E (■ a-c+d , b-c+d ; a+b’-G4d ; 1 -y ) . 

This giws Z 2 * substituted 

in (C-5), ve get 
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271 


OD 

/ 

0 


dz z“l'(2+o^, 


2+62; 25 -z) 



r(1 + 6'-a) 2 

( 2 f 

r(i+ ^ 2 '^ 


r (1 + a) 

r(i-a) 


a -♦ 1 


0 , 


(0-6) 



Appendix D 


SOLUTION OP THE APPROXIMATE INTEGRAL EQUATION POE 
SPECTRAL PUNCTION 


(a) Infrared region ((IWI- m) <: ) 

To solve Equation (2,16), we write it in the following 

form: 

(6)-1 )r ( w) ~5/d Cl? r( 0’ )(<)' ( w+u' ) 


. CO ( O - CO ' )' 

6' /d0*r(a)’) ^ T 

1 ’ri /m 


(D-la) 


a)^(co+1 )r(-co)£i 6 /dW r(-(o' )w' (co+co' )e(6o^-o«^) 


+ f 6' / do' r 

^ 1 • 


(-C0') 

p /m 


(L-lh) 


where t') > 0 and in the second terra on r.h.s. we have put 

w' 2 : 2, etc. Equations (D-la) and (D-1b) can be combined 

to write 


CO 


(cor^(6)) -^„(o)))- -6/da>' r )co' (to+co' ) e(o-2-(j ) 


.u2 


4 6' 5^ dco' r_(co') 'j " ■» (I)-2a) 

3-1 ti /m 


2 .2 


-id'? du'r ^ '/ • (Il-Zb) 

5 1 U /® ■ 
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Let us substitute Equation (2.18) in Equations (D-2a) 
and (D”2b): 

r' (w) -A rl(o)) + 2 6j^dw' Ar'(a)M 

1 

~ 6 /(jco' rf* (O' )(w-6» j (D-3a) 

Pf 1 

CO 

coAr^((j) - Ar^(co) + 2^/ <ico» Ar' (o') 

~ 4 6' ( ^' ) (o— o' )^, (L— 3b) 

1 + 

Next we substitute 

Ar^(o) = u^(o^) r^‘ (o) Ar^(o) = v^(o^) r°’(<'^) 

in Equations (D-3a) and (D-3b), differentiate twice and 
use the integral equations satisfied by r°* j then one can 
see that a possible solution is 

... 2 . ,..2, 6 ' r(“26) /,2 , ,2 . _ , 

U- ) - v. (“ ) = - -T "T ^ -1) + a (D-4) 

up to leading order in (o^-l). When substituted back in 
Equations (D-3a) or (D-3b), the constant of integration ’a’ 
turns out to be zero. 



(b) Asymptotic region 

let us perform the following break: 

r( 0 ) =e(w) (wr^Co)^) + r^(co^)). (p_ 5 ) 

Substituting Equation (D-5) in Equation (2.20) we have 



2 

m r 


-r 

c 

/co^) +. 

m (i'o((o^) 

~ (m^ ) ) ] 


n — 

Sj { 

-1 

^ w|w 

r-O) 

CO 

/ 

1 

) dm'co' 

e(m') (<i)r^(w'^)+ 

I'e(o)’^)) 


( / 

- Q 


m 



? 

(m+o)' ) (m'-wT 

+ 6' 


— 

/ , ) 

d(o’ r(w' 

)(m' (m+m' ) + 

P P 





6 m p /m 

"" *^2 

^.2 

i 

2 

dm’ 

(«r^(co'^) 

,2 

+ u' 





-6 ’ 

Ji 

m 

22 ^ 

(uir^(a^) 


(D-6) 


On equating odd and even parts in co, Equation (D-6) is 
equivalent to two equations. Again we break r(o) as 
r®(w) +Ar(w), where r°(m) is the solution of Equation (D->6) 
with li/m 0 and is given by Equation (2.21 ). Writing 

Ar^Cca) a rg(w^), A r^(co) = V 2 ( o? )r° (m^ ) 

and ecxpanding r® and r° asymptotically on both sides, we 
find that 
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U2(w )~ - 

O '(|j/io) U~^ + 

(P-Ta) 

Vg (Ci) ) 6 * 

(n/m)6 , 

,.2 2 

(B-7b ) 


Substituting back in the integral equations, we find b^=b2=0. 

To see how fax these asymptotic behaviors are applicable 
in the low-energy region, we proceed as follows. Consider a 
region such that 1 and H/m « (^-1) (one 

can think of p/m ~ 0l(>^-1) j). In this region we hare, 
from Equation (2.20), 


(^'^(Q- 1 ) r(w) +52/d0)' r (o)' )o)i (oh- Q ' ) 


Si -2 6» 


fA C<) — - 2 p 

( 7 - 7 ) (-] + I )r(w') (D-8a) 

11 ^ \i 


and 


r(-ti)) - dg -/dtij’r (-(i )o)' (c«><'W' 

w +|i/m Q P 2 

-2d’( / -OdQ' r(-a>*)(1+~2_ (tj-ax ) ), (D-8b) 

1 1 ■ ^ 


It is clear from Equations (D-8a) and (D-6b) that the 
contribution to r(+w) from the integral on r.h.s., is down 
by a factor of[ (y^m)/(w-1)] compared to the contri- 

bution from the l,h.s* The solution of Equations (I>-8a) and 
(D-8b) with r.h.s. equhl to zero is given by expression (2.21). 
Proceeding by the now familiar way, we write r(m) = r°(o)) +Ar(w) 
in Equation (Ma). Then in teims of even and odd functions 
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defined as in Equation (2,15»), {D-8a) reduces to 

-Ar_(o) + 26^ ^du’Ar„(co') 


Ci> -ix/xa 2 

~--2 6'( / - / )da)>(l + I )2 ^ 

11 ^ \x^ 

«+!.i/m to 2 „ 

-26' (/ " / .)dco' (1 + 4 ^^((0-0)' )'^)r°(-to» ) 

1 1 . 

2 -1-26 

-c26» / dco' (1 + j S- (co~<o' f) . (ylrl i... 

to-ti/m .. 2r(-262) 


— c 


22 


r(“262) 


1 

01 


I±lQ 


-1-2^2 


(D-9a) 


Similarly, frc»i Equation (D-eb), we get 



22 6 t X 

^ r (- 262 ) 


Jf 

m 


(m-1) 

T 


- 1 - 26 ^ 


(D-9b) 


With the substitutions 

o 2 o 

Ar_^(to) s Ar_,(co) = v^(o ) r^(to) 

in Equations (D“9a') and (I>*9b), one finds that a possible 
solution is 

Ug vl) « v^(o>^)k - ^ 6«(h/m) 


(D-10) 



Comparing Equations (]>7a) and (D-7b) with Equation (D-10), one 

2 2 

observes that one can identify ■u^(a> ) with U 2 (m ) in the 

2 

corresponding limit, but the same is not true for (co ) and 
2 

V 2 )• Thus it seems reasonable to write 

2 

A (6> ) u^ (m ) r°(m) X < |w} 00 . (2.23a) 

2 2 

On the other hand, (o) ) and Y 2 (c>) ) seem to be the limiting 

cases of some transcendental function. We shall write an 

2 0 

expression for A r_(m) which behaves as v^(w )r_((»)) and 
2 o 

V 2 (t^ ) r°(ci>) in the respective limits: 

A r_ (w )‘:i. - 5 ' (p/m) (w^-1 ) ^ur°((»>) ^ 00 , (2,23b) 



Appendix S 


APPROXIMATE INTEGEAi EQUATION POR SPSOTRAI 
PUNCTION ;.fflSN ^A/(W-m)«1 


We shall simplify expression for Im ^^[Squation (2.14b) J 

in the aboTe region neglecting terms of order 

l^W. We recall that I Wl , I ^ m . We also note that 

in the full integral Equation (2* "14), Im ^ ppears 

along with the spectral function P (W* ) whiob multiplies it, 

P (W ) is concentrated mainly at W’=«i. Thus terms like 2 » 

etc. appearing In the expansion of 1“ % oan be 

I'pp^'iiciLted by tfitb these approxi- 

-m (W' "m ) 

mations, simple algebra yields 


Im 


S(W,W')si--^ [ 6 WW' (W^-W'^) ©(W^-CW'+U) ) 
32 


+ 




When Equation (S -1 ) (along with Equation (P.l^a)) is substituted 
in Equation ( 2 . 14 ), in terms of dimensionless quantities, 
Equation ( 2 , 20 ) is obtained. 
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Appendix F 

SOLUTION OP THE GAP EQUATION 
USING ASIPIPTOTIC PREELOM 

Polloxo-Dg Lane's work^^, we shall solve the gap 
equation (3.10) using the fact that the QCD is asymptoti- 
cally free. 

We define a nonzero mass, in terms of which a 
XSB solution for the quark propagator may he written 

S“”’ (p,ii,tt^,^) P(p,p,mj^,gj^)-raj^ G(p,ii,ii^,gj^). (P-1) 

In terms of (P-1), the gap equation (3.10) reads 

G (p, p. ,11^ , gj^ ) (f ^ ) / (dk) ty ^»S (k,|-i,ii^, ^ ) ^ ^^ 

With K 9-s a dimensionless parameter used to scale 
the momentum p, the solutions to the renormalization- 
group (EG) equations for P and G turn out to be 

i’’(t<p,p,!%,gg) = exp(- P(p,ti,m(H), g(K)), (P-3a) 

and 

g ) 

G(Kp,ii,tflg,^) = exp(- I RTly Yja(x)]) G(p,p,m(H ),g(H )), 

% 


(i‘-3b) 
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whereas for the kernel the solution is 


exp (-2 i2a:||j)K(p ,ic,0;p,m(K ) ,g(K ) ), 

% 




where 


g^(K) ^ (t Ink-r"' + 


R-^oo 


(Inwr 


(5'-5) 


and 


m(H) 


K 


exp ( - / 

% 


g(K) dxYjjj(x) 


?T^ 


-) 


K > 




^■<,■*00 


(i’-e) 


For QCD with Nj, quarks, functions Yni^ ^ Landau 


gauge (for small g^^) are 

given 

by 


P Csh) = - I 4 > 

b = 

33-2 N™ 

(P-7a) 


c = 

1 

— 

(P-7b) 

Y (gj) ~ f 4 • 


\ 

(P-7c) 


Equation (P-2) can Tie rewritten in the form 




= (dk) 


(nkf • ♦ (Kky • • • Sj^ ^ 


(F-8) 



where 


K(p,k, ) -4 ^^*^5 ^cd ^ac ,d.b^^ ^ 5 ^ba* ^^”9) 

Erom Equations (E-3)-(E-g), it follows that G-(p,{a ,m(K ) ,g(K ) ) 
satisfies the same integral equation as G-(p,P ,Q^,gjj) ; 


G(p, l^,ci(V‘0jg(H )) 

G(k,ti,m(K),g(K)) K(p,k,0;p,m(ir<),g(4<)) 

= /(dk) —7--^ 5 5 . (5’-10) 

k E (k,. ..,g(l-0)-ni (h) G (,k^...,g(y)) 

We expect the m-even functions E and G to beha-ve in the 
limit ia(K) -«■ 0 as 

2 

E(p,^i,ra(K),g(K)) = E(p,|i,0,g(K)) + i in^(H) (^) 

am m=0 

+ 0(m^(K) In m^(K) ) (E-1la) 

G(p,|i,m(K), g(k)) = G(p,^i,0 ,g(k)) + 0(m^(K) Inm^OO). 

(E-"11b) 

A similar statement is true for the kernel as well 

K(p,k,0;lJ.,m(H),g(H)) = K(p,k,0; p,0,g(k )) +0(ra^(k)) . (E-12) 

Thus for H -00, Equation (E-IO) becomes 




G(kj P-i Ot , K{^ .. U K(p,k, 05 n, 0 ,g(K)) • 

k^E (k,p, 0 ,g(K)) 


(E- 13 ) 
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Taking advantage of asymptotic freedom, we approximate the 
kernel (in the landau gauge) by 


^ac ,db^^ ® ^ ^ — f^S (h) C2(5)(Y ) 


ac 


^(k-p);,(k-p)>'/(l!:-p) -Saj., , 

[ — : — yy:; )db 

(k-p ) + 


0(g^(K))j 


(i'-Ha) 


or, from Equation (P-9), 

K(p,k,05P,0,g(K)) = 3 iV(K) 02(3) 7"":2 ' + 0(g'^(‘<))- (5'-14b) 

(k-p) 


Our final approximation will be 


P^(k,p,0,g(K)) = 1 + 0(g^(K)) 


(P-IS) 


The integral equation (P-13) ultimately becomes (letting 
■p/\.i -* p and k/ \x -* k ) 

CKp,g{K)) = -3ig^K) CjCJ) / (dk) • <S-16) 

„ k (k-p) 

To solve Equation (P-16), we first rewrite it in euclidean 

form via a Wick rotation and then ope^^’’^© on both sides ty Q , 

Using the identities O — ^ — jr - -"4^ ^ (p-k) and 

(p-k) 2 

D G = 4(p^G) , where the derivative is with respect to p , 

Equation (P-1 6) can be written in the form 



99 


[ P 8(p,g(K))] j &(p,g(K))/p . (I'-n) 

1671 


The solutions of this equation are subject to the boundary 
conditions 


lim (p^ 
2 

p -► 00 


•' ""w ' + 1 ) G- = 0, 
dp 


lim 

2 

P ^ 



0 


G = 0. 

dp 


We define 


v= +11 -3g^(H) C2(3)/47 i^ 


(P-I8a) 


(P-ISb) 


fv» 1 - c/(b lnJ<), 


(P~19) 


where, from Equations (P-T), the ratio c/b is 


c _ 12 

F ” 33 - ^ Nj,' • 


(p-ao) 


The solutions to Equation (E-17), then, are 

c/b 

„ T) fv El ( "'l ^ 

+ 


4.1^ 

G+(P»g(H )) = a_^ p 


/InjL'i 

^+^ln*np'^ 


(E-21a) 


and. . 

(J„(p,g(K)) = a-p“”'“^2 “i (15^) 

P 

K '♦00 ^ 

Also G(Kp) is given, asymptotically, by Equation 


(E-21b) 

(E-3b) as 
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C/^Id 

G(Kp.,li,Tn^,gjj) ^[g^(H)J CG(p^0,g(^; )) + 0(m^ ( k) In m^(K )Ij, 

(F--22 ) 


Combining Equations (E-2la) and (E-22), we get 

■“c/b _2 g 

G^(hp,h, ^ % (InKp) + 0 (k (Inn) ). 

Here, B is an undetermined constant. Similarly, combining 
Equations (E-2lb) and (E-22), we get 

a , c/b o -R 

&-(KP,li,mR,g^) - ] + 0(H (InKr). (P-23) 

^ ^ “ P (Inn) 

Since this expression can be a function ef )f<; and p only 
through their product, must have the asymptotic form 
(up to a constant) 

a_ ~ (ln«)^^'^^ K (P-24a) 

(K ) ) exp [ -2/b g^ (k ) ] . (P-24b ) 

Accordingly, G_ will be written as 

S_Cp,h,i»b,Se) ~ +0(K’^(lnKp). (P-25) 

tt.=p/m[i - CD 

So far as we can tell from the BS equation alone, the error 
in G„ is of the same order as the solution itself. 

Firstly, we note from Equation (F— 22) that the 
asymptotic behavior G^(p)^(ln p)"^^^ exactly corresponds 
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to what we learn from a straight forward RG analysis, and 
does not differ from the case in which qiis-rks have nonzero 
bare mass and G(p) does not have purely dynamical origin. 

Thus we expect G_, which has a softer asymptotic behavior, 
to be the correct solution. This can happen if it corresponds 
to G arising solely from dynamics and is proportional, through 
the WT identity, to the bound -rstate vertex function 
P (P ,p+q) (at <1=0). To show this, we write 

ab 

iS(p)<p^(p,p+q) Y5 iS(p+q)6 

z 

= J' d^x (q)IT['i' (x) f (0) ]|0>, (P-26) 

^2 . - 

where trace over color space is understood. We may learn 

the asymptotic behavior of P , and hence of G, by taking the 
2 2 2 

limit Ip 1 -*■ 00, (p+q) /p -►I. Using operator-product 

expansion, one readily deduces 

9* Id 

iS(p)p (p,p+q)Y5 iS(p+q) 6 

~U(p,[i,mji,gg)<7r^(q)l ^ ( 0)1 Q>. (P- 27 ) 

Nonleading terms in (P-27) correspond to operators with 
canonical dimension four or more, and are expected to be 
down by approximately a factor of p relative to the exhibited 
leading term. 

Standard RG analysis of the Wilson coefficient function 


0 gives 
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U(Hp,[i,n^,gg) 



g(H) 

/ 

% 


dx 

PT^T 



(x)+y(x) ] } 


X U(p,4i,in(K), g(H)), (i'-28a) 

while from Equations (E-1), (E“3) and (P-6), we know that 

S-' (KP,., gj,) = Kexp[- ^)] 

% 

X s"'”' (p,p,m(k), g(/<)). (P-28'b) 


Finally from the WT identity of the renormalized axial 
vertex (in the cutoff theory) 

<l -2ni^(A.) ^ - 2^22 ”* (y^ S + S '* 

where is the renormalization constant appropriate tof*^ 
and m^ -*0 as /v oo , we see that m^ (A) 'Z,^/Z-q= 
and Zj^/Z2 are cut off -independent. This implies that 


It follows that 


(P*-29) 


<^(p,1i,niRTgR) oc- lim p (p,p+q) 
q -*0 


h=pAr 


cx> 


-2 ^ tor (x) 

, exp[ / T/i -T- j 


% 


WT 


(P-30) 

which is precisely the G_ solution. Moreover, we now know 

—3 g 

that the error in this solution is 0(H (InK) ). 
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Appendix G 


CHIRAI-SmiBlRY BREAKING SOLUTION K)E THE QUARK 
PROPAGATOR USING INPRARED PROPERTIES OP QCI) 


Here we shall derive results (3 .I 6) following 
Ref erence [ 35 ] . We start with the SB equation (3.13) 
satisfied by the quark propagator (3.14). The unrenormalized 
axial -gauge gluon propagator BjijXk) will be written as a 
sum of an infrared singular part given by (3.12) and an 
infrared finite part. For S(p-k)ry (p“k,p) S(p), longitu- 
dinal part given by Equation (3,15) will be used. With 
these substitutions, Equation (3.13) takes the form 


1 = 2 ^(|( P+G) - g^ 03 / (dk)YjiD,^a, (k) + 


p - P' 


- P’ 




+ Xp + |( terms. (U-l ) 

(T) 

In Equation (G-1 ) Xg, and are contributions of , 

Prom Equation (G-1) we can write equations for P and G 
separately which contain I and H respectively with (p.n) 
as coefficients. We ahAll specialize to the gauge p.n=0 
wherein H and I can be dropped and eiju^-tions for P and 
G take the sitnple forms 
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1 =p2E-g2o2/(ak)(J^6^„.-^[| 


p' p^+p;, p i) 


X D^^k) + }(j 


(G-2a) 


P 5’+g„C2Z(M)A ir/(dk) (1+ ) 

2(k.nr 


+ (F-F 


• 2 2 
-Tjit . .. 


2p.k-k 2(k.n)' 


+ Yj,, 


(G-2b) 


p^G-g^C9'^(dk) 


G - G 


P -P' 


Tt^|aa(p^“P.p' )+ P^P^ + PyP|^]D^y(k)+p^X, 


(G-3a) 


2 2 2 2 

= P^G+g^GpZ(M)A (dk) ^ (G-G')Cl + — '9 + 

° k"^ 2(k.nr 2 p.k-k" 


+ ] + Y. 

2(k.nr .® 


(G-3b) 


The functions Yp and Y^j axe and X^ plus the contributions 
due to nonsingular part of 2 ^ . 

lirst consider the following integral 

Ii 5 / ^ F ((p-k)^[ 1 + I (G-4) 

^ k^ ^ (k.n) - 


We shall assume that F is a well-defined function so that ' 

is finite. To deal with the singularity arising from k* 

we shall introduce a regulator prescription defined by the 
36 

replacement 

■* lim (lA^) (l/{k^+E^)). 

e -► 0 


(G-5) 



105 


Singularity arising due to the use of axial gauge will be 
dealt with, by using the principal value prescription^'^ 


— J ^ lim ^ [ Q 

(k.n; e» (k.n+ie*)^ 


+ 


Izll 




(-k,n+ie* 


J. 


(G-6) 


This prescription, in the present case, amounts to the 

replacBment 7sin© d9/cos^0 -*■ -2. With these prescriptions, 
0 

and after shifting the variable and angular integration, 
integral (G-4) becomes 

,,, 2 2 2,2 ..2 2 , 1/2 ,, 2 2 , 

T CD „2^,,2, +e ) -4k p ) ^ -Ik -p 1 ^ 

ii ^ f ) f- rr-5 7“ 

0 2p £ 2p 


TT^ 2 ^{2 .,£^2x1/2 ^ 

( (k +P ^ ) '"4k p ) 


(0-7) 


Here it can be observed that the order of taking the 
limit e -* 0 and perfcrming the integration is inportant. 
Thus if the limit is taken after the integration (for P=1 
in (0-4), ignoring ultraviolet behavior for the time being) 


T 2 2- r 

12=71 / dk L 

0 


{(k^+p^+^)^-4kV} ''/^-Ik^-p^I 


TTT 

2p e 


„ 1 

^ 

2p^ 


S 2 ' „ 2 2,1/2 ^ ^ 

{(k +p +e ) -4k p J 


whereas if the limit is taken before the integration 

.2 


2 

71- ^ j, 2 2 

‘■ w* / dk = — 7c , 

P 0 



This last result would have been obtained if the variable 
in (G- 4) were shifted without introducing e -regulator prescri- 
ption and the integration performed, We conclude that this 
difference is due to the presence of a 6 -function in the 
original integral (G-7) if the limit e 0 is taken before 


the integration, i.e., 



lim j- 
K-* 0 


P 


P ^ 9 o 9 9 9 9 

/ dk P(k'^) +TZ /2 2 dk “P ) 

0 P -K 

xp(k^)] . (G-8) 


That the r.h.s, of Equations (G— 4) and (G-8) are same can be 
checked by considering some simple algebraic functions for 
E, Thus if, for example, one chooses P((p-k) )=[(p“k) -p ]/ 
[(p-k)"^ + p^] so that the integral (G-4) is well-defined in 
the infrared region of integration (and prescription (G-5) 
is not needed), then the 6 -function part of (G-8) does not 
contribute indeed as it should. 

However, relation (G-8) can not be applied when 
E = constant. The reason is that in this case the original 
integral (G-4) is logarithmically divergent in the radial 
variable while the angular integration of the bracketed term 
is zero (Equation (G-8) takes only the latter into account). 

To do this kind of integral we shall use dimensional regulari- 
zation in addition to the two prescriptions mentioned earlier: 
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-r* _ , r. 1 

II - / [1 + 2- 


2 2 


k n 


(k.n)‘ 


,2<*^ 

"* liin 7 ~ " . H". '. " M ^ 

0 lc^(k^+e2) 


[1 


+ 


1 i2 2 

1 k n 


Second term can be written as 


64 


(k.n) 


2" 1 >2 Cl) < 4. (G-9) 


n 


PV / 


d k 


• (l?T?)(k,n)‘ 


2 o )-2 


2t»), 
d k 


k (k +e ) 


(G-10) 


Prom Equations (G-9) and (G-10), we get 
4-2« , d^^ 

I = X 

k (k +6 ) 


^ ^ci) g2(a».2) |^P(2-a)) 

2 ■' 

> % • 

lira (jy* 2 

e •* 0 


(G-1 1 ) 


Obviously non-zero contribution (G-11) has come from the 
ultraviolet region of integration. 

Next, consider the following integral 


o o 2 , 2 

^ ^P((p-k)^)[^£-^ 




2p.k-k 


r? (2p. Arfe^.). 
2 (k.n) 


]. 


(G-12) 


Again with the usual assumption that P is a well defined 
function and employing the regulator prescriptions mentioned 
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earlier, after shifting the variable and performing the 
angular integration (G-12) takes the foim 


J = ^ dk^ P(k^)[ - — (k^+p^+e^)^-4kS^}*''''-lk^-p^l 

0 2e"^(p -k ) 


1/2 


1 


— 2 

P -k 4p 


(4 .^) .^J^l'(p2,l,^e2)2-4kV 

4p4p-k4 


1 2 2 2 , 2 . 
k +p . p -k r 1 


4p'"(p^“k^) 2£?^ fe^-p^l {(k^+p^+e^)^-4k^p'^' 

(G-13) 


]L 


Prom our experience of dealing with the integral (G~4), we 
expect that the first derivative of 6 -function, in addition to 
6-function itself, will occur when the limit e -* 0 is 

taken inside the integral. No higher derivative of 6 -function 
should occur. This can be seen by taking an explicit form of 
P, like P((p-k)^) = ((p-k)^-^^)V((p-k)^+p^)“, n=2,3,...., 
in which case no contribution should come from 6 -function or 
any of its derivatives. Thus if the limit e -*■ 0 is taken 
inside the integral (G-13), then it will take the form 


lim 271? [ 

K 0 


2 2 

0 


P(k^)^+/ dk^P(k^) 
.2 2 

P -K 


2p 


{ a 6' (k^-p^)+b6(k^-p^)}] 


2 2 . 


(G-14) 


where the derivative of 6 -function is with respect to k , 
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Explicit ca-lculation with a couple of simple algebraic functions 

P Q 

for P(k ) in Equations (G-I 5 ) and (G-14) shows that a = 2p'^ 
and b = -1 . Thus, finally (G-12) reduces to 


P(k^) -271^ P(p^)-4/p^ 


(G- 15 ) 


Here again we can not apply Equation (G-I 4 ) for P = constant. 
However integrands of (G-12) and (G-4) have the same ultra- 
violet behavior except for the difference of overall sign. 

Hence their ultraviolet contributions will cancel each other 
in Equations (G-2b) and (G-3b) for P = J'(p^) and G(p^) 
respectively, and hence Equations (G-8) and (G-14) are sufficient 
there for P = constant, as well. Using Equatior^ (G-8) and 
(G- 14 ) in Equations (G-2b) and (G-3b), we got 


1 = p P(p ) 4- pM p 


^2 mvf) 


+ ) 


(G-16) 


0=pW)+|^ J _fif-G(p2) 


+ pM^p^ ^ 




(G-17) 


Where p=AG2Z(M) 


Analogous equations for the renor- 


-r 

malized quark propagator can be written with the alteration 


!■ ' S - Sg, Tj - Yj,g (1-Yj, = Zjd-YjgWand Yg - Y^g, 
Next, we shall drop and and denote renormalized 

functions without the subscript R, Then the equation for P 
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can be conTerted into the following second order differential 

2 2 

equation (y = -p /^M ): 

j P" + (1-y) pt - P = 0 . (G-18) 

Independent, solutions of Equation (G— 18) are the confluent 
hyper geometric functions §( 1 ; 1 ; -p^/pM^) and S ' (1 ; 1 j-pVpM^). 
When substituted back in Equation (G- 16) (with Yj, = 0 ), first 
solution is ruled out while for the second one, the multi- 
plicative constant is fixed: 

P(p2) = 1^- W ( 1 ; 1 j -p^/^M^). (G-19) 

PM 

Similarly, the equation for G can be converted into the 
following second order differential equation: 

y G" - (y - 3/2) G» - 2G = 0 . (G-20) 

Independent solutions in this case are $ (2 » 3/2 ; -p^/p ) and 
( -p^/pl^ )'"'''^^ § ( 3 / 2 - 1/25 -p^/pM^), When substituted in 
Equation (3-14), second solution gives, 

effective mass = S ^(0) = G '^( 0 )'= 0 , 

and hence is ruled out. Thus we get 

G(p^) = ^ (2;3/2; -p^/pM^). 


(G-21) 
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Appendix H 

SOME ANGULAR INTEGRATIONS USED IN CALCULATION 
OR PAULI POEM PACT OR 


We shall perform the following integration: 


(y^^-2. )y [6 

(k-p)^ k'^-k A p. TJ- 0 


A0 


n. (k-p) 


(k-p)x(k-p)cyn 

^ J ^ 

(n.(k-p) ) 


(H-1) 


Pirst consider 


d Q. 


I, = / 


k 


. (k-p)'^((l‘^+21c.q) 




(k-p)xnari- (k-p) on \ 
n. (k-p) 


(k-p)x(k-p)0n 

+ n J. 

(n. (k-p)) 


(H-2) 


Choose the gauge n,p = n.q = 0 and assume that p q , 

r H' 

then 


(k-p)^ (q^+2k.q) 


k r . 271^ 

^ L 2r , 2""2- 

k +p +pq k>.Ik -p I 


1 ^ + g - Ik ^ l l } j(h-3 ) 


+ K — i 

k +p^+pq k > 


2 ' 2 2 

where we have used (k> = max (k ,p )) 
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d Q, 


(k-T))' 


27^ r ^ ^k 

k^> ’ TZZT^ 


(k-p) 


271^ 

k> Ik^-p^l 


d 0, 
k 

“2 

q +2k.q 


d Q 


P/ 


k 


q +2k.q 


= 4[1-(1 - 4 -)^/^ e(#-k2)] 

k q ^ 

= S(k,q). 


In the last integral (which is not defined for 2k >q) we 
have taken the principal value because in the context of 
the whole integral (H-1), it is the principal value which 
is relevant. 

Por those integrals in which (n,k) or its power 
occurs in denominator, we shall use the standard principal 

Ca 

value prescription which in simple cases reads 


'1 1 


(k.n) ' ^ { [(k.n) +ie'] ^ +[(k.n) -ie* ] 


-1 


='0 


(H-4a) 


(k.n) 


(k.n) + ie* +[(k.n)~ie' } 


n Sin© d© ^ 

f = -2, 

o Cos © 


(H-4b) 


Thus, 


d Q, 

“k 

(k-p)^(q^+2k,q)(n.k) 


= 0 


kx 


(k-p)^(q^+2k,q)(n.k) 


T 

T 

n 


(H-5) 
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and 


n 

(k-p)^(q +2k.q)(n.k)^ 


- 471 ^ k 

n y L— n « — ' - n T? + — n <■ —n 

k +p +pq k I (k “p )^| k +p +pq k Jk^-p | 


1 -I 

r k + -p ,1 ; 1 

L ' "M - — -1- "- VS ' "?s t " ■■ n"" -— n" ' 


+ 


9(9^4 - f ) j ] 

2 ' 4 2 2 v 1/2 

Kk (q -4k q ) '^ • 


(H-6) 


where we have used 


dG, n 

f k ^ 

(k-p)^(n.k) 


.^2,, 2 2^ 
A-Ti (k +p ) 

Ic Kk -p ) 1 


dOj ^ n ^ 

(k-p) (n.k)^ 


A ^ 

4 It 

k Ik “p I 


dOj^ n 

— 2 2 

(q +2k.q)(n.k) 


dG n 

p / — g ■ • ' ' n - 

- (q +2k.q)(n,k) 


47 t ^ e ( qV 4-- k ^) 

, 2 , 4 " T 2 ^ Ti /2 • 
k (q^-4k q ) ^ 


In the last integral we have again taken the principal value 
for the same reason as explained above. 


(k-p)'^(q''^+2k.q) (n.k)^ 


27i2(k2+p2)q^ I ^ q ^ 

n ^ q ^ k^l ( k ^- p ^)^ 2 n ^ 


(H-T) 


d 9 ky^ k^j 


F X — rt 15 = A q-xq-+ B nxn-+ 0 6^ 

(k-p)^(q‘^+2k.q)(n.k) 


(H-8). 
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A — ^ ^ r —1 — ^ -r , -z ^3 •] 

2q. n n p 

1 I. ,2 I 

2n n n p 

1 I-t v2 

G=-5'[“^“ “2‘^2'^“^^' 
n n p 


(H-8a) 


^,,2 2 . ,, 2 2 x 2 

[i - g.( k +p , ) + ..(jsltplL.] 

- (q‘'+2k.q)(n.k) (k-p) (k-p)"^ ' 


^ j. ^ ^k 


_ 7^^Q(q^/4-k^) ^ pW . _4 


'^Q(q^/4-k^ 


■7" 


{ 


1 

^ — 5 — T 


2n k +p +pq k Ik -p I 


,, , } ■ ii±vb l . 


lo. 


.k2(,^-4kV)'/^ - 4n2 

(H-8b) 


Equations (H-3), (H-5)“(H-8b) are substituted in 
tion (H-2), we get 



<l^‘a 

a 


1 


Mj . ■ e(p2-k2) 

2p2(k2+p2+p,)2 


+ ’‘pH(| . p . ^+?p<i) e(k%^) 

2p (k +p +pq) 
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71 


7 p / 4 2„2.1/2 

^ — 5 — 2 2" ^ ^ — i — ^ — T 

(k +p +pq) q^(k +p +pq)^ 


+ 


— 4 A A~^ ^2 — 2 — 2 J^(<a^/4-k^) ]. 

(q^-4k'^q^)^/2(k%p%pq) • 


(H-2') 


Next, we consider 


d Q, k 

(k-p)'^(q^+2k.q) 


k [2 + a.- {k: . £j- 3 

(n.kr - 




(H-9) 


where 


A. 


2/[- 6 ,( 2^ - ^ I _3L 9(p^-k^) 

k p 471 : 2p'^(k +p +pq) 


2 2 , 2 , , „ _ 2 , 4 2 2 j/2 

. i^^ H(8r: . ±4p . q l 9(k%2) ^ ^ { ( . q.-4k p ---p 

47 r 2p (k +p +pq) q (k +p 4-pq) 


+ 




2 2 — : ^ e(qV4-k^) I 

^-4k q ) (k +p +pq) - 


(H-9a) 


is not of our interest. Brom Equations (H-2), (H-2'), 
(H-9) and (H-9a), we get 
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^ (Y >(-2k )Y [6, - (^-V>a’^X 

(k-p) k' -k n. (k-p) 


(k-p)x(k-p)j3n , 
^ J 

(n. (k-p)) 


f ai,2 ”^>''i8PgHgq) g; 

0 2p (k +p +pq)' 


00 p 

-■*■ / dk 
2 
P 


2fiL(8E +4p.^CTXk^j„ - 271^ ^ dk^{ .La. ^- .^ ^ ..... ^ . , ) - , ^ .. g , fe ... 

2p (k +p +pq) ^ <1 (k +p +pq) 


2G(k^) 


(qMk'"q^)^/‘"(k^+p^+pq) 


}]- 27t^ ^ [ / dk^ G(k^) 


f 1+4^ (8p^+4TDq) 

p ■ 4n: 2p (k 4-p +pq)^ 


^ ? dk^ 

r\ A ^ 


( 8 p^+ 4 pq) G(k^) 
2p (k +p +pq) 


2 


.4,y7v2a2J/2^y^2. 


+ / dk t “^■*5 — — ^ .. v- i. 

0 q (k +p 4-pq) 


+ y- ■ T > 

(q - 4 k q ) ^ (k +p +pq) - . P 


(H- 10 ) 


Next, consider the following integration 


(k-p) 


G(k«2) Y V OY w r;^ (k-p)y,n^-f (k-p),jn^ 

^ ,2 _^^2 Y ^"2k ji)Y ^ ^ - n. (k-p ) 


(kT))x(k-p)^yn 

■f J 

(n. (k-p)) 



= / 


+ 


where 



Here, 






(k-p' )xn^(k-p' )an 


n.k 


(k-p' )x(k“p' )cy n 

-2 J +o -function part (if any), (H-11 ) 


(n.k) 


p' = p+q. Also define p^^ = w' ci|_i. Then, 


dQ, 


r -o- 

_ « / j 


(k-p' ) (q +2k.(-q) ) 


[6 


(kvp' );>,ncy+(k-p' )afnx 


X(S 


n.k 


(k-p' )\(k-p* )dn 
+ ] 


(n.k) 


= - ^xu^P**P'j ^P •* “<1P*, K •* -K‘ } 


^ -4 [ ■ e(p'2^c2) 

q k 2p''^(k +p' -p'q) 


.2 


_ 2^2 ,2 , .2 p ^ ^ TT^ 

2p' (k +p' -p'q) (k +p’ -p'q) 


, 4 ,,2 2 . 1/2 2 o o 

- 2 ( ■ y + TT 7 T - i/2 ' "2 2~' »(4 


q (k +p' -p'q) (q -4k q )' (k +p''^-p'q) 


(H-12) 


2 2 2 2 

: +p' -p'q = k +p +pq. It can be easily checked that 


dl.2 4^2 ^0 . 


Hence by shifting the variable, as in Equation (H-11), the 
apparent (infrared) divergent nature of integral does not 


introduce any 6 -function. 


In a similar way, 


M- 


d Q V k,, , .2 2 q n., 

= ^ ^ [ 2 + A* -^ + B' 

(k-p' )(2k.q-q ) (n.k) - ^ q'^ ^ n 


(H--13) 


where 


A:j = A,|(p-*p', qp-*-qp», h- “H' ) 


k 4^ 2p» (k +p^ -p^q) 


^ ^ ”/’j8p'"-4p . 'a ) e(k2-p,2) + a! 


"'-W 

4lt 2P’ (r+p' -p'qr 


,4_.^2„2ll/2 


^ .. y,^ ^ ; 1 

q (k +p’ -p*q) (q -4k q )'' (k +p' -p’q) ■ 


X e(qV4-k^)]. 


(H-13a) 


Substituting Equations (H-12 )-(H-13a) in Equation (H-11), 
we get 


d\ 


(k-p)^ k' 




n.-(k-p) 

2 

(k-p)x(k-p)an ^ 

(n. (k-p ) ■ 


4 - 
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. P' 


0 2p' (k +p +pq)^ p’ 


" 5(1^") 


4- 


2 ,1 ^ li ^ P ,, 2 „ /, 2 . 


7TTr^‘"27T727r2~'2 — ■/ [■ / dk G(k ) 
(q -4k q ) '^ (k +p +pq) - 0 


x ■2 li>^: . (8p'^"4p*q} > 

2p'‘='(k^+p'^+pqr 


? dk^GCt^) i i ^'Lsf^^P't) 
P’ 2p' (k +p + pq) 


q^/4 


2% T 
0 


dk^{G(k^) I j 

q (k +p. ^-pq)"^ 


2G(k^) , 

+ y "■2r7'"’2"2\i75-,; 2’"' '2 — 

(q -4k q ) ^ (k +p +pq) - 


S 2 

— 

q p' 




D * 2 2 

^ dk G(k ) + n terms, 
0 


(H-14) 


Substituting Equations (H-lo) and (H-14) in Equation (H-1), 
we get 


(k-p)^ 



^x(Y^,J^-2k^)Y^[6^^ - 


(k-p)\(k-p)on^ 

(n. (k-p))^ - 


(k-p)^ng+ (k-p)^^ 

n,(k-p) 


«■» tVy ,Y„JtA' /' ak‘ 


S(k^ )j2p^^-p 


X, ‘ 

P' (k 


7 

+p +pq) 


^ -k' 


2 

q A 2 2 

/ dk G(k )x 
0 


q (k +p +pq) (q‘^-4k^q^ (k^ +p^+pq ) ' 
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/ dk^G(k^) ^ J dk^ G(k^)]+ n teims. 

0 P 0 - 1^ 


(H-15) 
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